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Abstract 

We present a unified framework for the preclusion of non-degenerate multiple steady 
states in a network of interacting species. Interaction networks are modeled via systems of 
ordinary differential equations in which the form of the species rate function is restricted by 
the reactions of the network and how the species influence each reaction. We characterize 
the set of interaction networks for which any choice of associated rate function is injective 
within each stoichiometric class and thus cannot exhibit multistationarity. Our criteria 
rely on the determinant of the Jacobian of the species rate functions that belong to the 
class of so-called general mass-action kinetics. The criteria are computationally tractable 
and easily implemented. Our approach embraces and extends much previous work on 
multistationarity, such as work in relation to chemical reaction networks with dynamics 
defined by mass-action or non-catalytic kinetics, and also work based on the graphical 
analysis of the interaction graph associated to the system. 

Keywords: influence specification; reaction network; monotone kinetics; Jacobian; degener- 
ate; general mass-action; interaction graph 

1 Introduction 

Interaction networks are used in many areas of science to represent the structural form of a 
dynamical system. This is in particular the case in systems biology and biochemistry where 

biochemical reactions arc represented in the form of a network. However, similar network 
structures are also used in ecology, cell biology and epidemics, as well as outside the natural 
sciences, to describe the possible interactions or reactions between some species of interest. 
Common to these interaction networks is that they consist of a set of species and a set of 
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interactions among the species. The state of the system is given by the concentration (or 
abundance) of each species and each interaction represents a transformation of the state of the 
system. An example is the chemical reaction A + B ^ 2C where one molecule of A and one 
molecule of B form two molecules of C. 

Typically, a system of ordinary differential equations (ODEs) is used to describe how species 
concentrations change over time. The rate function of the system describes the instantaneous 
change in the concentrations when considering simultaneously the individual rates of all re- 
actions in the network. Reaction rates are generally unknown but some qualitative aspects 
might be assumed and inferred. For instance, the presence of the species on the left side of a 
reaction (A, B in the above example) is a prerequisite for the reaction to take place and higher 
concentrations of these species typically lead to higher reaction rates. In some cases, reaction 
rates are fixed to follow a specific functional form that might depend on parameters to be 
inferred from experimental observations. In other cases, only weak assumptions are imposed 
on the reaction rates. We consider restrictions given by a so-called influence specification,!^ 
which specifies how each species affects (positively, negatively, or neutrally) the reaction rates. 
In this work an interaction network refers to the network (species and reactions) together with 
an influence specification. 

For many networks the structure of the interactions and the influence specification alone 
determine dynamical and steady-state properties of the system (for example, multistationarity, 
persistence, or oscillations). That is to say, irrespectively the rate functions and the parameters 
quantifying them, taken together with the initial species concentrations, the system shows 
qualitatively the same type of behavior. It is perhaps surprising as the network structure itself 
does not encode any information about the specific rate functions and abundances. Even small 
networks might have many parameters which potentially could give rise to a rich and varied 
dynamics, as well as differences in the long-term behavior of the system. 

Of particular interest has been to determine whether a system allows for multiple positive 
steady states, also known as multistationarity. Multistationarity provides a mechanism for 
switching (rapidly) between different responses and confers robustness to the steady-state values 
of the system.H^'iSl Pqj- example, the growth (response) of a bacterial system depends on the 
availability of nutrients (stimulus). If the system exhibits multistationarity then fluctuations 
in the supply of nutrients determine the growth behavior. Within the range of nutrients where 
two (stable) steady states exist, the steady-state values are robust to minor fluctuations in the 
stimulus. 

One way to address whether a system exhibits multistationarity is by finding the positive so- 
lutions to the steady-state equations of the system. Solving the equations might prove difficult, 
if not impossible, with difficulty depending on the assumptions about the reactions rates and 
the number of parameters. Here we take a more conceptual route and focus on understanding 
the characteristics of interaction networks that cannot exhibit multistationarity, irrespectively 
of the specific choices of reaction rates. 

Various criteria have been developed to preclude the existence of multiple positive steady 
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states for general classes of reaction functions, also called kinetics, such as mass-action kinetics,^"^'^'^ 
non-catalytic kinetics,'^ and weakly monotonic kinetics These criteria typically utilize the 
structure of the system together with some assumptions about the form of the rate functions. 
For example, mass-action kinetics assume that the reactions take place in a well-stirred con- 
tainer such that the frequency of collisions between species is proportional to their abundances. 
The reaction 2 A + B ^ D would have a reaction rate proportional to the product of the species 
concentrations, c^c^, and the steady-state equations become a system of polynomial equations. 
Capitalizing on the polynomial form of the equations has lead to specific conditions to preclude 
multistationarity.'^ ^ ^ 

The aim of this paper is to provide a unified framework for injectivity of a network while 
allowing the kinetics to be of a general form but respecting an influence specification. Injec- 
tivity refers to injectivity of the (species formation) rate function that governs the dynamical 
behavior of the system. If this function is injective for the allowed kinetics then the system, 
or network, does not have the capacity for multiple positive steady states. The idea of in- 
jectivity was introduced in Ref. [3| but it is also underlying previous work on the preclusion 
of multistationarity.'^ We show that injectivity of an interaction network is closely related to 
injectivity of an interaction network taken with general mass-action kinetics, introduced in 
Ref. [llj. General mass-action kinetics are a generalization of mass-action kinetics and con- 
fer greater flexibility in the form of the rate functions than mass-action kinetics. Ref. [21] 
emphasizes the importance of general mass-action kinetics and the closely related fractal ki- 
netics in biochemistry but their importance outside biology for modeling purposes is also well 
documented.lf^ Both are based on a power-law formalism and it can be argued that they 
provide approximations to kinetics in general.'^ The latter is also exemplified in our work. 

General mass-action kinetics have many things in common with mass-action kinetics and 
parallel results can be derived for the two types of kinetics. We derive necessary and sufficient 
determinant criteria for a network to be injective over different classes of general mass-action 
kinetics, defined by influence specifications. The determinant refers to the determinant of a 
modified version of the species formation rate function (Definition 5.2 in this paper). We show 
that the determinant of the Jacobian of the modified function is non-zero for all concentration 
vectors and kinetics in one of the classes if and only if the network is injective over the particular 
class. The results utilize and extend parallel results for mass-action kinetics."^ The determinant 
criteria depend in part on the network structure and in part on the influence specification. 

These criteria are of importance for several reasons. First of all, they provide means to pre- 
clude multiple positive steady states for classes of general mass- action kinetics. Secondly, they 
are easily computable using symbolic software packages such as Mathematica or Matlab and, 
thus, they are of practical use. Importantly, the set of interaction networks that are injective 
for all choices of general mass-action kinetics agrees with the set of interaction networks that 
are injective for all choices of kinetics that are in accordance with the influence specification. 
Therefore, our determinant criteria determine injectivity more broadly. 

This work extends and embraces previous determinant criteria developed for networks taken 
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with mass-action kineticsP'^'^ In Ref. [23], the authors also characterize networks that are 
injective for classes of kinetics defined by an influence specification (defined in a slightly different 
way then in this paper). The criteria discussed in this paper are different from the criteria in 
Ref. [23] and we extend their results. Further, our work clarifies the role played by general 
mass-action kinetics in deciding injectivity. 

In Ref. [2| the authors study injectivity of a broad class of kinetics by putting some restric- 
tions on the form of the reactions. Namely, reactions of the form 5 — )• belong to the network 
for all the species S and kinetics are the so-called non-catalytic kinetics. Interestingly, under 
their assumptions one of our determinant criteria for injectivity follows from their criterion. 
Further, in Ref. |2] injectivity of the network is derived by providing a criterion that guarantees 
that the Jacobian of the species formation rate function is a P- matrix (defined in Section 10). 
We extend this result to other classes of kinetics and arbitrary networks and show that the Ja- 
cobian of the modified species formation rate function is a P-matrix if and only if the function 
is injective. 

Finally, this work also relates to a criterion for multistationarity based on the interaction 
graph.nSl'211'phe interaction graph records the sign of the entries in the Jacobian of a dynamical 
system and does not assume the existence of an underlying network structure as in our case. 
We show that their criterion agrees with one of our determinant criteria applied to a certain 
type of interaction networks, namely, those consisting exclusively of reactions of the form 
5 — )• and — )• S". We show that the existence of multiple non-degenerate steady states in the 
class of dynamical systems considered in Ref. \16\ [M] is precluded by our injectivity criterion 
for these interaction networks. It follows that the results presented here are stronger since 
the information provided by the network structure enables to preclude multistationarity for a 
broader class of dynamical systems. 

The structure of the paper is the following. In Section 2 we introduce some notation and 
in Section 3 we introduce the basic concepts of networks and kinetics. In Section 4, we discuss 
interaction networks and give examples from the literature. Section 5 introduces the notion of 
degeneracy and injectivity, two key concepts. Section 6 focusses on general mass-action kinetics. 
In Section 7, we derive a number of results on injectivity for different classes of general mass- 
actions kinetics. These results are used in Sections 8 and 9 and extended to broader and more 
general classes of kinetics. Section 9 includes also the relationship of our work to Ref. [23] . In 
Section 10, we develop the relationship to the work of Ref. [2j and in Section 11, we relate our 
work to that of Ref. [16j . Finally, in Section 12 we show that other types of kinetics could be 
used in place of general mass-action kinetics. To keep the exposition clear in the main text, all 
proofs are in the Appendix. 
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2 Notation 



Let M+ denote the set of positive real numbers (without zero) and M+ the set of non-negative 

real numbers (with zero). Similarly, let N be the set of non-negative integers. Given a finite 

g 

set <5, let N be the semi-ring of formal sums v = J2e&£ ^eE, with A^; e N. If A^; € N for all 

E £ £, then we write v S N^. The semi-rings and are defined analogously. 

The ring of polynomials in £ is denoted The total degree of a monomial ITsef E^^ , 

with ue a non-negative integer for all E, is the sum of the degrees of the variables, Y2e££ ^e- 
The degree of a polynomial is the maximum of the total degrees of its monomials. 

If a polynomial p vanishes for all assignments a: £ ^ M_(_ then p = identically. Further, 
if p is a non-zero polynomial in Wy£] such that the degree of each variable in each monomial 
is either 1 or zero, then all the coefficients of p are non-negative if and only if p[a{£)) > for 
any assignment a: £ ^ M-|_ . 

For a vector v = {vi, . . . , Vm) £ W^, let the positive part be v+ = (max(?;i, 0), . . . , max(r;m, 0)) 
and the negative part be V- = (min(wi, 0), . . . , min(fm, 0)). The support of v is defined as the 
set of indices for which v is non-zero, supp(?;) = {i\vi ^ 0}. The positive support of v is 
supp"'"(?;) = supp(u+) and the negative support is supp~(u) = supp(t;_). Let denote the 
transpose of v and u ■ v the usual scalar product in M". 

For every x € M, we let sign(x) G { — 1,0, 1} be defined as 



sign(x) 



ifx<0, 

if X = 0, 

1 if X > 0. 



3 Networks as dynamical systems 

In this section we introduce networks and kinetics, and associate a dynamical system with a 
network and a kinetics. The definition of a network is identical to that of a chemical reac- 
tion network, which is used mainly in (bio)chemistry to describe networks of (bio)chemical 
reactions.'^ In general we use the nomenclature that is standard for chemical reaction networks. 
See for instance Ref. [6] f7] for background and extended discussions. 

Definition 3.1. A network J\f consists of three finite sets: 

(1) A set 5 of species. 

(2) A set C C N of complexes. 

(3) A set 7^ C C X C of reactions, such that (y,y) ^ TZ for all y £ C, and if y G C, then there 
exists y' £ C such that either (y, y') £ TZ or (y' , y) G IZ. 
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A network is denoted hy M = {S,C,Tl). 

We use the convention that an element r = (y, y') GTZis denoted by r : y ^ y' . The reactant 
and the product (complexes) of a reaction r: y ^ y' are y and y' , respectively. By definition, 
any complex is either the reactant or the product of some reaction. The zero complex G C is 
allowed by definition. A reaction 5 — )• 0, S E 5, is called an outflow reaction. 

Throughout the paper, we use n to denote the number of species in S. We fix an order 
in S so that S = {Si, . . . , Sn} and identify N with N". The species Si is identified with the 
i-th canonical n-tuple of N with 1 in the i-th position and zeroes elsewhere. Accordingly, a 
complex y £ C is given as y = Y17=i Vi^i iuii ■ ■ ■ ^ Vn)- We assume that y E N as reactions 
typically involve integer numbers of species. However, the results presented in this paper hold 
generally for y E M .In examples we will often use other letters than Si for species to ease the 
presentation. 

Example 3.2. Enzyme biology provides a rich source of examples. For instance, consider the 
network with set of biochemical species S = {S*!, 5*2, Ss, 5*4, ^5, ^e}, set of complexes {Si + 
53, Si + 5*4, 5*2 + 5*4, 5*2 + 5*3, 5*5, 56 } and reactions 

Si + 5*3 < ' 5*5 ^ Si + 5*4 5*2 + 5*4 < ' Sq ^ S2 + 53. 

That is, the reactions are ri : <S'i+S'3 — )• Ss, r2 : Ss — )• 51+5*3, r-^: S^ ^ 5i+54, r4 : 52 +54 — )• S%, 
r^: Sq ^ S2 + S4 and rg : 56 — )■ 52 + 53 . This network is a main building block in protein 
modification systems and is known as the futile cycle. It assumes the Michaelis-Menten enzyme 
mechanism in which a substrate 53 is modified into a substrate 54 through the formation of an 
intermediate 55. The reaction is catalyzed by an enzyme 5i. The modification can be reversed 
via a similar set of reactions with an intermediate 56 and an enzyme 52 . 

Reactions in a network are schematic representations of dynamical processes. Over time 
the concentrations or abundances of the species in the network change as a consequence of the 
reactions. In order to describe the dynamical properties of the network we introduce a kinetics 



(Definition 3.3) and the species formation rate function (Definition 3.7). The kinetics provides 
the reaction rate for given species concentrations and the species formation rate function the 
instantaneous change in the concentrations when considering simultaneously the rate of all 
reactions. 

Definition 3.3. A kinetics for a network J\f = {S,C,TZ) is an assignment to each reaction 
y — )• y' E 7^ of a rate function Ky^yi : Qk ^+ , where Qk is a set such that M" C ^ 
and 

Ky^yi{c) > for all c E VIr- 

A kinetics for a network is denoted hy K = {Ky^y/)y^y/^fi and the set of kinetics for J\f 
by 1C{N). If Ky^yi is differentiable for all y — t- y' E 7^ and c E M" then K is said to be a 
differentiable kinetics. The set of differentiable kinetics for Af is denoted by ICd{M). 
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Example 3.4. Let M = (5, C, IZ) be the network with S = {Si, S2, S3}, C = {Si + 5*2, 53} and 
TZ given by the reaction Si + S2 ^ S3. The kinetics K £ K.{M) defined by Ks-^+S2^S3{c) = 
kci/{{f3 + 01)02), where k,a,/3 are positive constants has ^Ik = M x M-i- x R-|-. 

Example 3.5. Kinetics commonly used in chemistry and biology are the so-called mass- action 
kinetics. These were introduced by Guldberg and Waage in the 19th century based on the ideal 
assumption that the rate of a reaction is proportional to the product of the concentrations of the 
reactant species P Specifically, each reaction y — )• y' is assigned a positive constant ky^y' G M_|_ 
and the rate function for the reaction is given by 

n 

Ky^y' (C) = ky^yl Cf 

i=l 

with Q.K = 1^+ and K = {Ky_^yi)y_^yi^Ti. Under in vivo conditions, however, the use of mass- 
action kinetics might not be fully justified. Reactant species might not form a homogeneous 
mixture, for instance because they appear in too low concentrations or because their distri- 
bution depends on spatial constraints. In situations in which the use of mass-action is not 
justified, the use of other types of kinetics such as general mass-action kinetics or Hill-type ki- 
netics are often preferred. These will be introduced later. Note that Ky^yi{c) is an increasing 
function in q and does not depend on other species than those involved in y. 

Definition 3.6. The stoichiometric subspace of a network M = (5, C, TZ) is the following vector 
subspace of M": 

r = (y - y| y y G TZ). 

Two vectors c, 0' £ M" are called stoichiometrically compatible if c — c' G T, or equivalently, if 
uj ■ c = OJ ■ c' for all oj G T-*" . 

We let s be the dimension of T and, thus, the dimension of is d = n — s. Note that F is 
independent of the choice of kinetics and depends only on the structure of the reactions. Being 
stoichiometrically compatible is an equivalence relation which partitions M, into classes, called 

— n 

stoichiometric classes. In particular, the stoichiometric class of a concentration vector c G M_|_ 
is {c + r} n M^. Knowing the value of w • c for all vectors a; in a basis of T-*- determines the 
stoichiometric class to which c belongs to. The value of cj • c is called a conserved amount. 

Definition 3.7. The species formation rate function for a network M = (S, C, TZ) with kinetics 
K G K.{M) and stoichiometric subspace F is the map fx '■ — ^ F defined by 

/x(c)= Ky^y'ic){y'-y)- 
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The dynamics of a network N with kinetics K G 1C{N) is described by a set of ordinary 
differential equations (ODEs) given by the species formation rate function: 

c = /ii-(c), 

where c = c{t) denotes the derivative of c{t) with respect to t. The steady states of the network 
are the solutions to the system of equations in ci , . . . , obtained by setting the derivatives of 
the concentrations to zero: 



=/x(c) 



This system of equations is referred to as the steady-state equations. In particular, we are 
interested in the positive steady states, that is, the solutions c to the steady-state equations 
such that all concentrations are positive, c G . Observe that the image of Jk is contained 
in r and it follows that for any u) S T"*- we have that • c = 0. Therefore, w • c is independent 
of time and determined by the initial concentrations of the system. 

For any differentiable function / = (/i, . . . , fq): — )■ M'^ defined on a set O including Mip, 
let Jc(/) denote the Jacobian of / at c G M^, that is, the q x m matrix with entry (i, j) being 



dfi{c)/dcj. If K e K.d{J^) is a differentiable kinetics then the Jacobian Jcifx) '■ 
evaluated at c € WL is given for any 7 = (71, . . . , 7„) E Mi by 



roil 



K 



E E 

i=l y^y'eTZ 



dci 



{y -y)- 



(3. 



IS 



Example 3.9. The stoichiometric subspace of the futile cycle in Example 3.2 

r = (5*1 + 53 — 5*5, 5*5 — S*! — 5*4, 5*2 -I- 5*4 — 5*6) 

and has dimension s = 3. The dimension of F-*- \s d = n — s = ?> and a basis is {w^, w^, w'^} 
with 



.1 



51 + S5, u:' = S2 + Sq, = 53 + 54 + ^5 + 56. (3.10) 

If X G fC{M) is any kinetics then the corresponding system of ODEs for the futile cycle is: 

c = Ks,+ss^ss (c)(-l, 0, -1, 0, 1, 0) + Ks,^s,+s,{c){l, 0, 1, 0, -1, 0) 
+ Ks,^s^+sAc){l, 0, 0, 1, -1, 0) + Ks,+s,^s,.{c){0, -1,0, -1, 0, 1) 
+ Ks,^S2+s,{cK0, 1, 0, 1, 0, -1) + Ks,^S2+s^{c){0, 1, 1, 0, 0, -1), 
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which written component-wise yields the system 



Cl = -Ks^+Ss^SrAc) + Ksr^^Si+sAc) + Ks^^si+Si{c) 
C2 = -Ks2+S4^Seic) + Kse^Si+sAc) + Ksf^^Si+Saic) 
C3 = -Ks.+Sa^Ssic) + Ks„^Si+sAc) + Ks,,^S2+S:i(.c) 
C4 = Ks^^^Si+sAc) - Ks^+Si^Seic) + Kse^Si+sAc) 
C5 = 7^:51+53^55(0) - i^55^5,+53(c) - Ks^-^Si+sAc) 
C6 = Ks^+s^^Saic) - Ksa^Si+Siic) - Kse^S2+S3{c)- 

Observe that ci + C5 = C2 + ce = C3 + C4 + C5 + ce = for any kinetics K G K.{J\f). 

4 Interaction networks and influence specifications 

Here we introduce interaction networks and some concepts related to them. An interaction 
network consists of two parts: a network and an influence specification. The latter stipulates 
structurally what species influence the rate of a reaction and whether a species acts as an 
inducer or an inhibitor of a reaction. We use influence specifications to restrict the set of 
possible kinetics for a network. Our definition of an influence specification is inspired by 
Ref. [23]. 

Definition 4.1. An influence specification for a network M = {S,C,TZ) is an assignment to 
each reaction y — t- y' G 7^ of a function ly^y' '■ 5—7- { — 1, 0, 1} . An infiuence specification for J\f 
is denoted by X = {Iy^y')y^y'e'iz- We define the following distinguished influence speciflcations: 

• The complex dependent influence specification, denoted by 2^:, is defined as the infiuence 
specification with Iy^yi{Si) = 1 for z E supp(7/) and Iy^yi{Si) = for i supp(y), for all 
y^y' ^n. 

• The reaction dependent influence specification, denoted by T-ji, is defined as the influence 
specification with Iy^y/{Si) = 1 for i E supp(y), Iy^yi{Si) = —1 for i E supp(?/') \supp(y) 
and Iy^yi{Si) = for i supp(y') U supp(?/), for all y ^ y' £ TZ. 

The species Si with Iy^yi{Si) / influence the reaction y — y' . Those with Iy^yi{Si) = 1 
are assumed to have positive influence or enhance the reaction while those with Iy^yi{Si) = — 1 
are assumed to have negative influence and an inhibitory effect on the reaction. Those with 
Iy^yi{Si) = have neutral influence. 

Remark 4.2. In typical applications an influence specification satisfies 

ly^y'iSi) i^Q, for ah iEsupp(y). 
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Negative influence 



Positive influence 




Figure 1: Reaction dependent influence specification for the futile cycle, drawn as a bipartite 
graph. There is an edge between a species and a reaction if the species has non-zero influence 
on the reaction. If the influence is positive, the edge is solid. If the influence is negative, we 
draw a dashed edge. The complex dependent influence speciflcation is in the case obtained by 
removing the dashed edges. The species interacting in a reaction, for example S5 — )• S*! + ^4 
cannot be read off from the edges. 



That is, the enhancers and inhibitors include all species involved in the reactant complex y. We 
do not impose this restriction here because the condition is not required for the mathematical 
development of the work. In chemical reaction theory it is typically required that the species in 
the reactant complex have positive influence on the reaction?^ I^Sl and not negative or neutral, 
as in Definition 4.1 However, relaxation of this assumption is found '^1 ED ^nd if some species 



in some reactions are only partially known, as might be the case in some biochemical systems, 
this appears reasonable. 

Definition 4.3. An interaction network (J\f,I) is a network J\f = {S,C,TZ) together with an 
infiuence specification I. 

If no species appear both in the reactant and the product of a reaction then the reaction 
dependent infiuence specification Iti is given by the functions ay^yi{Si) = sign(yj — y^), y — >■ 
y' G TZ. However, if this is not the case then Cy^y', y ^ y' £ TZ, might not define the 
reaction dependent infiuence specification. For instance, the reaction Si + S2 ^ Si + S3 gives 
ay^yi{Si) = 0, contradicting that IniSi) = 1. 

It is sometimes useful to illustrate an infiuence specification with a labeled bipartite graph 
with node set 5 U 7^. We draw a positive edge between a species and a reaction if the species 
has positive infiuence over the reaction. We draw a negative edge if the species has negative in- 
fiuence on the reaction. An example using the futile cycle and the reaction dependent influence 
speciflcation is illustrated in Figure [T} 
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G2 

Figure 2: The influence specification for the transcription of the three genes are shown. The 
presence of G3 reduces the production of Gi , whereas the presence of Gi and G2 cooperatively 
induce the production of G3. Likewise Gi induces the production of G2- 

Example 4.4. Common examples in the literature are gene regulatory networks.'^ These are 
typically represented by diagrams as the one shown in Figure [2j The diagram represents three 
genes, each transcribing a protein. The proteins mutually affect the transcription rates of the 
genes such that the genes influence each other. The diagram corresponds to the interaction 
network with production reactions 

^Gi i = 1,2,3, 

and influence specification X = {Iy^yi)y^yii=.']i with (only non-zero values are listed) 

h^Gi{Gz) = -1, and /o-s>G2(G'i) = Iq^gAGi) = h^Gz{G2) = 1- 

In the following, we assume that S, C, and TZ denote the sets of species, complexes and 
reactions, respectively, of a network AA, unless otherwise specified. Let {M,X) be an interaction 
network. For y ^ y' £ TZ, define 

Iy^y> = {i\Iy^y'iSi) = l}, ly^y, = {i\Iy^y>{S^) = -I}, and I^^y, = {i\Iy^yiSi) = 0}. 

Two concentration vectors a,b £ are said to be non-overlapping with respect to I, 

if Iy^y> ^ supp(a) imphes C supp(6) 

for all reactions y — )• y' . That is, the coordinates ai,bj (potentially with i = j) cannot both 
be zero if i,j are both in Iy__^yi for some y — )• y'. When it is clear from the context what 
infiuence specification we are referring to, we omit "with respect to X" and just say that a, b 
are non-overlapping. Non-overlapping is a concept that specifies how two concentration vectors 
on the boundary of the positive orthant are positioned with respect to each other in relation 
to an influence specification X. In particular, if one or both of the concentration vectors a, b 
are positive then they are non-overlapping. 
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Definition 4.5. A kinetics K = {Ky^yi)y^yi^'ji for a network M = {S,C, TZ) respects an 
influence specification X if, for all c G ilx, 

Ky^yi{c) > if and only if ly^y, C supp(c). 

If a kinetics respects an influence specification then a reaction occurs if and only if all species 
with positive influence on the reaction are present (that is, they are in positive concentrations). 
If one species is not present then the reaction cannot occur. Absence of species with negative 
or neutral influence does not prevent the reaction from taking place. 

Let ^K{Iy^yi) ■= {c G ^k\ ly^yf ^ supp(c)} = {c S R \ Ky^yi{c) / O} dcuote the set of 
concentration vectors for which the kinetics Ky^yi does not vanish. 

Definition 4.6. Let K = (Ky^y')y^y'^-ji be a kinetics for a network M = (S,C,TZ) that 
respects an influence specification I. We say that K is strictly monotonic with respect to I 
if for all reactions y ^ y' £ TZ and i = 1, . . . ,n, the restriction of the function Ky^yi{-) to 



(i) strictly increasing in Cj if i G 



(ii) strictly decreasing in Cj if i G 



(iii) constant in Cj if i G ly^yi ■ 

Let /C(A/', X) denote the set of kinetics that are strictly monotonic with respect to the influence 
specification X. 

The definition says that the rate functions Ky^yi are strictly monotonic or constant in 
the coordinate Cj whenever the remaining coordinates take positive values for the species with 
positive influence. Strictly monotonic refers to the kinetics, but this does not imply that the 
species formation rate function is monotone. 

Example 4.7. The length of a larvae is often assumed to increase linearly with a slow down 
in the growth rate as the length increases. Denoting by c the length of the larvae, one model 
of the length is c = ai(a2 + c), where ai,a2 are positive constants. This system arises from 
the interaction network with reaction 5* — t- 25 (the addition of new cells at the two ends of the 
larvae) and kinetics Ks^2s{c) = ai/(a2 + c). This kinetics is strictly monotonic with respect 



to the influence specification given by Is-^2s{S) = —1 (compare Remark 4.2) 



Lemma 4.8. Let K = {Ky^yi)y^yi^-ji he a kinetics for a network M that respects an influence 
specification X. Then, K is strictly monotonic with respect to X if and only if for each pair 
of non- overlapping concentration vectors a,b G Qk, the following implications hold for all 
reactions y ^ y' : 



12 



(i) if Ky^yi{a) > Ky^yi{b) tkeTi Si & S with sign(aj — 6j) = Iy^y'{Si) 7^ for some i. 

(a) if Ky^y,{a) = Ky^yi{h) thcu either Ui = bi for all i G I^^y, U ly^yi, or sign(aj - 6j) = 
Iy^yi{Si) 7^ and sign(aj — bj) = —Iy^y'{Sj) ^ for some distinct i,j. 

Example 4.9. A standard model of population growth i^^ 

c = rc(l-^y (4.10) 

where c denotes the size of a population S and r, K are positive constants. The system can be 
considered an interaction network M with two opposing reactions S ^ 2S (birth) and 5" — t- 
(death), S = {S}, C = {5, 25", 0}, and influence specification Is^2s{S) = Is^o{S) = 1. The 
kinetics with Ks^2s{c) = rc and Ks^o{c) = r(? jK is strictly monotonic with respect to X and 



gives equation (4.10). The system has a stable steady state at c = and an unstable steady 
state at c = 0. The two steady states are non-overlapping with respect to X. 

This population model has been generalized in various ways, also incorporating competition 
among different populations Any model of one population where c initially increases in c until 
a certain point and then decreases can be seen as giving a strictly monotonic kinetics for M . 



One particular instance of a competing species model is shown in Example 4.13 



Definition 4.11. A kinetics K = {Ky^y/)y^yii^n for a network A/" is dijjerentiable with respect 
to an influence specification I if K respects I and, for every reaction y ^ y' G TZ, Ky^y'{-) is 

continuous at c G i^K, differentiable at c G M", and for each i the partial derivative -—§^^{c) 
has constant sign 



W(5.) = sign(^^(c)) 



in M". Let /Cd(A/',X) denote the set of kinetics that are differentiable with respect to the 
influence specification I. 

Note that /Cd(AA,X) is not the intersection of K.{J\f ,X) with ICd{M) as we require the sign of 
the partial derivatives to be constant in Definition |4.11[ which is not implied by being strictly 
monotonic and differentiable. Any mass-action kinetics K belongs to /Cd(AA, Xc). We state 
without proof: 

Lemma 4.12. If K £ K.d[M ,X) then K £ 1C{M,I). 

Example 4.13. An example from ecology is given by the Lotka-Volterra equations for modeling 
predator-prey dynamics or competing populations.'^'^ The modeling equations are 

ci = ci(a - /3c2), and C2 = -£2(7 - /3ci), 
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where ci and C2 are the abundance of two species Si (prey) and 52 (predator) . The system can 
be considered a network M = {S,C, TZ) with S = {81,82}, C = {0, 81, 82,281,282, 81 + 82} 
and n = {8i^ 28i, 81 + 82^ 282, ^2 ^ 0}. Define a kinetics K £ }C{Af) by 

Ksi-^2Si{ci,C2) = aci, K 3^^+82-^232(^1,02) = /3ciC2, and Ks2~^o{ci, 02) = 7C2, 

with a,/3, 7 > 0. Then, with the complex dependent influence specification Xc, {N ,Xc) is an 
interaction network and K G /C(A/',2c) is a mass-action kinetics (Example |3.5[ ). The species 
formation rate function is given by fxici, C2) = ( ci(a — /3c2), —02(7 — /3ci) ). The stoichiometric 
subspace 

r = {81, 82 — 81, —82) 

has dimension 2. Hence T"*" = {0} and there are no conserved amounts. The system has one 
stable steady state (ci, C2) = (7//3, a/ (3) and one unstable (0, 0). It shows periodic oscillations.l^Sl 
The reactions with the given kinetics can be interpreted in the following way (here assuming 
prey are mice and predators are owls): the mice reproduce, the owls reproduce while feeding 
on mice (one mouse is consumed for each newborn owl) and owls die. Mice only die when 
consumed by owls. The mouse population finds ample food at all times and reproduce at 
a rate proportional to their abundance, whereas the reproduction rate of the owls depends 
linearly on the availability of mice. Lotka and Volterra proposed the model independently of 
each other as a model of species interaction.'^SI In epidemics, the same model has been proposed 
as a model of disease dynamics with prey being individuals susceptible to disease and predators 
being infected individuals .^^^^ Also the model has appeared in physical chemistry as a model 
of H2O2 combustion.!^ 

Remark 4.14. The number and form of the reactions, as well as the influence specification, 
is not uniquely determined from the equations describing the system. In Example |4.13 the 



reaction 5i — )• 28i could be replaced by — )• /Si, but the interpretation of the reaction might 
change. Similarly, the second reaction 51 + 52 — ?■ 252 could be replaced by 52 — )• 252, assuming 
that species 5i has positive influence on the reaction. 



Example 4.15. Ref. [13J used a model similar to that of Example 4.13 to study bacterial 
uptake of nutrients in microbial ecology. The modeling equations are 

V^CiC2 , . VciC2 

c\ = — , and C2 



K + C2' Y{K + C2)' 

where Y, K, V are positive parameters and ci , C2 the concentrations of two species 5i , 52 where 
5i is a bacteria that feeds on a (chemical) nutrient 52. The system can be represented as a 
network A/" = {S,C,1Z) with S = {5i,52}, C = {0, 5i,52} and IZ given by 

5i *-25i, 52 0, 

representing the clonal reproduction of the bacteria and the consumption of the nutrient. The 
kinetics K with Ks-^^2S-i (ci, C2) = and Ks2^q{ci, C2) = y|x+c2) ^® strictly monotonic with 
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respect to the influence specification X with ls^-^2S\{S\) = Isi-^2Si{S2) = 1 and Is2->o{Si) = 
Is2^o{S2) = 1- All steady states have ci = or C2 = 0. Hence two steady states are non- 
overlapping with respect to I if at least one of them is strictly positive. Further, all steady 
states are degenerate. 



5 Degenerate steady states and injectivity 

In this section some key concepts and definitions are introduced, namely that of degeneracy of 
a steady state and injectivity of a network. 

Definition 5.1. Let M = {S,C,TZ) be a network and K G ICd^M) a differentiable kinetics. A 
steady state c G M" of A/" if degenerate if ker:{Jc{fK)) n F / {0}. 

Alternatively, a steady state is degenerate if the Jacobian restricted to the stoichiometric 
subspace F is non-singular. Denote the components of the species formation rate function by 
Ik = ifx,!, ■ ■ ■ , fK,n)- 

Definition 5.2. Let M = {S,C,1Z) be a network. A basis {uj^ , . . . , u)'^} of F-*- with cj* = 
(A^, . . . , A^) is said to be reduced if A^ = 1 for all i and A*- = for all j = 1, . . . , i, . . . , d. If 

K € 1C{N) is a kinetics for N then the associated extended rate function fx '■ — ^ I^" is the 
function defined by 

/x(c) = (W^ -C, ...,tj'^-C,/x,d+i(c),...,/A',n(c)). 

After reordering of the species set if necessary, a reduced basis always exists for any network 
M and is unique. Therefore, from now on, we fix that the set of species is ordered such that a 
reduced basis exists. It does not depend on the kinetics K £ K,{J\f). 



Example 5.3. The basis of F provided for the futile cycle in equation (3.10) is reduced. The 



associated extended rate function fx'- — ^ ^ for any kinetics K £ 1C{N) is 

/a-(c) = (ci + C5, C2 + C6, C3 + C4 + C5 + C6, Ks.^^Si+Si{c) - K s^+Si-^Seic) + Kse^s-i+sAc), 

The next proposition is proved in Ref. [8] for mass-action kinetics. However, the proof 
applies generally to differentiable kinetics without changes. 

Proposition 5.4. Let N = {S,C,TZ) be a network, K G ICd{M) a differentiable kinetics, 
{u!^, . . . ,uj'^} a reduced basis ofT^ and fxic) the associated extended rate function. If c £ M" 
then 

ker('/c(/K))nF = {0} if and only if det(Jc(7i^)) / 0. 
In particular, a steady state c G M" is degenerate if and only if det^Jcifx)) = 0. 
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Finally, we introduce the notion of injectivity. 
Definition 5.5. Let N = {S,C,TZ) be a network. 

(i) We say that M is injective over /Co ^ 1C{N) if for any pair of distinct stoichiometrically 
compatible vectors a,b £ M" we have fKio) 7^ fxib) for all K £ ICq. 

(ii) If X is an influence specification for N, we say that N is I-injective over /Cq ^ IC{M,I) 
if for any K £ ICq and any pair of distinct non-overlapping stoichiometrically compatible 
vectors a,b G Qk we have /a- (a) 7^ fxib)- 

A network is said to have the capacity for multiple positive steady states over /Co if there 
exists a kinetics K £ ICq and distinct stoichiometrically compatible vectors a,b £ M" such that 
fx (a) = fK{b) = 0. Clearly, if M is injective over some appropriate set of kinetics /Cq, then Af 
does not have the capacity for multiple steady states over /Cq- If an influence I for M is specified 
and M is X-injective, then the existence of pairs of distinct non-overlapping stoichiometrically 
compatible steady states a,b £ Q,k is also precluded. Therefore, the concept of non-overlapping 
vectors accounts for the occurrence of some types of multiple steady states on the boundary of 
Oft- and not only multiple steady states in the interior of Qk- 

Note that in Definition |5.5K ii) above, the concentration vectors a, b are required to be non- 
overlapping with respect to I. This set of non-overlapping pairs might vary with K since i^K 
depends on K £ fCQ. Since pairs of positive concentration vectors are non-overlapping for any 
influence specification I, X-injectivity implies injectivity. In Ref. [23] , the authors preclude 
the occurrence of pairs of distinct stoichiometrically compatible steady states a, b such that 
at least one of them is in the interior of (which is taken to be M_|_). As noticed above, 
such pairs are non-overlapping and hence covered by our approach. In Ref. [8], the condition 
of non-overlapping is applied to networks with mass-action kinetics and influence specification 
Ic- It is possible to have fx {a) = fxib) for a pair of not non-overlapping vectors and at the 
same time fxio) 7^ fK{b) for all non-overlapping pairs. ^ 

The rest of the paper is devoted to characterize the networks that are X-injective for different 
families of kinetics and that consequently cannot have the capacity for multiple positive steady 
states. The aim is to provide a criterion for a network N to be X-injective over a set of kinetics 
^0 ^ /C(AA,X) in terms of computational tractable quantities that rely exclusively on the 
network structure and the influence specification. To this end we introduce the class of general 
mass-action kinetics (defined in Section [6]) and derive some injectivity results for classes of 
general mass-action kinetics using techniques introduced in Ref. [8j. Furthermore, in Section [sj 
we will relate injectivity of a network with general mass-action kinetics to X- injectivity for the 
same network with kinetics that are strictly monotonic with respect to X. 
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6 General mass-action kinetics 



General mass- action kinetics form a large family of kinetics ."J^ They are generalizations of mass- 
action kinetics and are based on a power-law formalism. Their general form make them flexible 
for modeling purposes in many areas of science such as chemistry, ecology and epidemics. As 
mentioned in Example |3 . 5 1 general mass-action kinetics might be used when mass-action kinetics 
are not suitable and, furthermore, they can be seen as an approximation to a wide range of 
other kinetics;-- Furthermore, these kinetics are important in the study of injectivity in that 
they, in some sense, are "dense" in the set of kinetics that are strictly monotonic with respect 
to an influence specification (to be made precise in Section [s]). That is to say, injectivity of a 
network M over certain sets of kinetics can be determined on the basis of injectivity of Af over 
suitable sets of general mass-action kinetics. 

In this section we introduce networks with general mass-action kinetics. Networks with 
general mass-action kinetics have properties that are very similar to networks with mass-action 
kinetics and similar results regarding injectivity can be derived for general mass-action kinetics 
as for mass-action kinetics themselves. 

For a concentration vector c and v G M", we associate the power-law = ]Xi=i c^'- We 
define O" = 1 and 0^ = for any real number x ^ 0. For example, ii v = (2.1, 0.72, 0, —1) G M^, 
then the associated power-law is c" = 02'^'^ c^^ . A power-law is well defined for all c S M_|_ 
such that Q > if < 0. 

Definition 6.1. Let N = {S,C,TZ) be a network. We define: 

(i) A kinetic order is a 7^-tuple of vectors v = (vy^y')y^y'i=-Ti with Vy-^yi G M". 

(ii) A rate vector is a 7^-tuple of scalars k = {ky^y')y^y',=n with ky^yi G M-|-. The constant 
ky^yi is called the rate constant of reaction y ^ y' . 

A kinetics K = {Ky^yi)y^yi^n G 1C{N) is a general mass-action kinetics with rate constant k 
and kinetic order v if 

Ky^y>{c) = ky^yic^y^y' , fov all y y' G 7^, 

and is the maximal common domain of definition of the power-laws c"y^y' , y ^ y' £ TZ 
in the positive orthant. For simplicity, we identify the pair (k, v) with the kinetics K. We let 
/Cg(AA)[v] denote the set of general mass-action kinetics {k, v) with arbitrary k but fixed v and 
let ICg{M) denote the set of all general mass-action kinetics for the network M. 

Note that ky^yi is a scalar while Vy — ^yf IS On vector. By definition, general mass-action 
kinetics are differentiable kinetics, that is ICg{M) C jCdi-M). Mass-action kinetics are special 
types of general mass-action kinetics obtained by considering the kinetic order with Vy^yi = y 
for all y ^ y' G TZ. Therefore, if we let y = {y)y^y>^ii then the set of mass-action kinetics for 
a network M is K.g {M) [y] . 
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If K = (k, v) G lCg{M) is a general mass-action kinetics for a network M then the species 
formation rate function is denoted by fx = fK,v and, for c G ^k, we have 



The Jacobian of /^^v with respect to c G M" takes the form 

Jc{fK,^r){l) = X] ^j/-*-j/'c''^^^'K->2/' *c7)(y'-y), (6.2) 

y-^y'en 

with u *c = UiWi/ci. 



Remark 6.3. If k G is not fixed then the function fK,v{c) can be seen as a polynomial 
function in the variables ky-^yi for all y — > y' G 7?.. 

A kinetic order for a general mass-action kinetics is intimately related to an influence 
specification for the network. If v is a kinetic order for J\f, define an influence specification 



Iy-^y'{Si) — < 



1 if i G SVlPV'^ {Vy-^yl) 

— 1 if i G SVlY>^~ {Vy^yl) 

otherwise. 



Reciprocally, if X is an influence specification for A/^, let e(X) = {ey^yi{X))y^yi^Tz be the kinetic 
order defined by 

^y^y'i^) = i^y-^y'i^l)^ ■ ■ ■ ^ ^y-^y' i^n)) ■ 

Note that X(e(X)) = X. If I = X(v) then we say that I is the influence specification associated 
to the kinetic order v. 

For an influence specification I, let K,g{M ,X) C K,{J\f ,X) be the set of general mass-action 
kinetics that are strictly monotonic with respect to X. If v is a kinetic order then it is straight- 
forward to see that any kinetics («, v) G Kg{M) is differentiable with respect to the asso- 
ciated influence specification X(v). Therefore, ICg(M,X) C lC,j(J\f.X). Likewise any kinetics 
(k, v) G K,g{M) that is strictly monotonic with respect to an influence specification X fulfills 
X = X(v) . Hence, we have the following lemma: 

Lemma 6.4. For any kinetic order v we have /Cg(A/')[v] C ICg{M,X{v)). A kinetics (k,v) 
belongs to Kg (A/", X) if and only if I(v) = X. Further 

ICgiJ\f,X)= U ICgiMM. 

v|X(v)=X 
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Example 6.5. Consider the futile cycle of Example 3.2 The kinetic order y = (yi)i=i,...,6 of 
any mass-action kinetics is {yi denoting the vector corresponding to reaction r,): 

yi = (1,0,1,0,0,0), y2 = (0,0,0,0,1,0), = (0, 0, 0, 0, 1, 0), 

y4 = (0,1,0,1,0,0), ys = (0,0,0,0,0,1), ye = (0, 0, 0, 0, 0, 1). 

Assume, for instance, that the concentration of the demodification enzyme 52 acts as an en- 
hancer or inhibitor of the modification of S3 by the enzyme Si, that is, of reaction ri. This 
effect can be accounted for by considering a kinetic order v = (ui)i=i,...,6 of the form 

= (l,i;,l,0,0,0), ^;GM\{0}, Vi = i > 1. 

In this case, the species formation rate function /k,v is given by (denoting by ki the rate 
constant of reaction r^): 

Ksi+Ss^s^ic) = kiCiC3C2 Ksr^^Si+Ssic) = k2C5 Ks^^Si+Siic) = ksc^ 
^52+54^56 (c) = /C4C2C4 Ks(i^S2+Siic) = k^CQ Ksa->S2+Sa{c) = kQCQ. 

If for instance v = 0.5, then Ksi+s^^Ssic) = kiciC3C2'^ while iiv = —0.5, then Ksi+Sa^Ssic) = 
kiCiC3C2^'^ . The influence specification associated to v, X(v), has the following non-zero signs 
(Ij refers to reaction rj): 

/i(Si) = l /i(52) =sign(^;) Ii(53) = 1 12(^5) = 1 13(35) = l 

h{S2) = l 14(^4) = 1 15(^6) = 1 16(^6) = 1. 

The only possibly negative influence is /i(5'2). Note that X(v) depends only on the sign of v. 
In this example, the kinetic order e(X) is given by 

ei = (1, sign(^;), 1, 0, 0, 0), 62 = (0, 0, 0, 0, 1, 0), eg = (0, 0, 0, 0, 1, 0), 

64 = (0,1,0,1,0,0), 65 = (0,0,0,0,0,1), 66 = (0,0,0,0,0,1). 

Example 6.6. So-called SIR (S=susceptible, I=infected, R=recovered) models are standard 
in epidemiology to describe the outbreak of an epidemics in a population. One particular SIR 
model^ considers the network with set of species {S, I, R} and reactions 

S + I ^21, I ^ R. 

The first reaction says that a susceptible individual might become infected in the presence of 

an infected. The second reaction says that infected individuals eventually recover. The SIR 

dynamics can be expressed in different ways. One possibility is the following set of differential 
equation^211 

ci = — A;ic^C2, C2 = kiCiC2 — ^262, C3 = A;2C2, 
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where ci, C2, C3 are the concentrations of the species S, I, R, respectively, and /ci, A;2 > 0, / 
are the parameters of the model. That is, Ks+i^2i{c) = kic'(c2, and Kj^ji{c) = /C2C2. The 
amount N = S + I + R is conserved. All steady states have C2 = and C3 = N — ci, which 
depends on the initial concentration of S. Hence, all infected individuals eventually recover. 
Further there are always multiple positive steady states in any stoichiometric class (given by 
N). 

The associated influence specification X is given by 
Is+i^2i{I) = Ii^dl) = 1, Is+i^2i{S) = sign(z.), Is+i^2iiR) = Ii^dS) = Ii^r{R) = 0. 

The parameter 1/ accounts for inhomogeneity in population mixing. If z/ = 1 then the population 
is homogenous and the disease spreads according to the law of mass-action, whereas if 7^ 1 
then the population is inhomogenous and the kinetics is a general mass- action kinetics. If is 
negative then susceptible individuals are repelled by infected individuals. If we allow = 0, 
then the spread of the disease would be independent of the number of infected individuals. 

Many variants of the SIR model have been proposed,'^ not only at the level of the ODEs 
but also at the level of the interaction network. For example, individuals could be allowed to 
die, or recovered individuals could become susceptible again. 

7 Injectivity for networks taken with general mass-action ki- 
netics 

The characterization of injectivity for networks with general mass-action kinetics will be key to 
provide a characterization of X-injectivity for networks with general kinetics. In this section we 
provide criteria for a network to be injective with respect to different (nested) sets of general 
mass- action kinetics. 

7.1 Fixed kinetic order 

We first study kinetics with a given fixed kinetic order. 

Theorem 7.1. Let N = {S,C, TZ) be a network and v a kinetic order. Then the following are 
equivalent: 

(i) M is injective over JCg{M)[v]. 

(ii) ker(Jc(/K,v)) n F = {0} for all c G M*^ and k M^. 

(Hi) det(Jc(7«,,v)) / for all c(^Wi^ and k e M^. 

To simplify the notation we introduce the following: for any set of m reactions, R = {yi — )■ 
y'l, ■..,ym^ y'ni] ^ 7^, let 
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• Z{R, v) be the n x m matrix whose i-th column is Vy.^yi 

• r(i?) be the n x m matrix whose i-th column is y'^ — yi 

• For a set J C {1, . . . , n}, let Z(R, v) j and T{R)j be equal the matrices Z{R, v) and T{R), 
respectively, with the rows j £ J removed. 

The next proposition is an adaptation of Ref. ^81 Prop. 6.3, Cor. 7.3]. Let A/" be a network 
and let 0{M) = {i| — )• IZ} be the set of indices for which outflow reactions do not belong 
to TZ. Finally, let Od{J^) be the set of subsets of 0{M) of cardinality d and let 1 = (1, ... , 1) 
denote the vector with 1 in each entry. 

Proposition 7.2. Let M = {S,C, TZ) be a network, v a kinetic order, and s the dimension 
of the stoichiometric subspace. The determinant dei{Jc{fK,v)) is a homogeneous polynomial in 
K = {ky^yi)y^y'^Ti of total dcgrcc s and linear in each rate constant ky^yi . 

Further, let R = {yi y[, . . . ,ys ^ y'^} be a set of s interactions from TZ. The coefficient 
of the monomial ^^=1 ^Vi^y' det( Jc(/k,v)) for c £ M" is 

^-i+EU-,,^,^ ^ det(Z(i2,v)j)det(r(i?)j)nc,- 

Remark 7.3. The determinant is in general not a polynomial in v or in c because the coordi- 
nates of V enter the expressions as exponents of c. However, the determinant det{Z{R,v)j) is 
a polynomial function in the non-zero coordinates of Vy^yi, y ^ y' £ R, excluding the entries 
with indices in J. 



Remark 7.4. Note that the expression in Proposition 7.2 differs by a sign (—1)* compared to 



the expression in Ref. [8, Cor. 7.3] where the coefficients are given in terms of —T{R)j. 
Proposition 7.5. Let J\f = {S,C,TZ) be a network. The following are equivalent: 
(i) M is infective over K.g{J\f)[v\. 

(a) The non-zero products det(Z(i?, v) j) det(r(i?) j) have the same sign for all sets R of s 
reactions from TZ and J G OdiM). Further det(Z(i?, v) j) det(r(ii)j) ^ for at least one 
choice of R and J . 

The matrix Z{R, v) j only relates to the influence specification X on A/" and not the structure 
of reactions. Similarly, the matrix T{R)j only relates to the structure of the reactions. 

Example 7.6. Consider the futile example with the kinetic order v introduced in Example 
16.51 We have 

det(Jc(/K,v)) = -cicl{c2 + vc^ + C4)kik^ki - ciC2kik^k^ - cici^kikska 

- cl'^^{ci + C3)kik4kQ - 02/02^4^6 - C2/C3^4^6- 
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If f > 0, then all the terms of the determinant expansion have the same sign and by Theorem 7.1 
the futile cycle is injective over ICg{M)[v]. If, on the contrary, v < 0, then the term t = 
vcid^c^kik^kii has sign opposite to the rest of the terms. It follows that the futile cycle is not 
injective over /Cg(A/')[v] if t> < 0. The term t corresponds to the set of reactions R = {^i, rs, 
and J = {1, 3, 6}. Indeed, for these sets we have 



Z(i?,v)j= 1, T{R) 





so that det(Z(i2, v)j) = -v, dei{T{R)j) = 1 and c-^+^«=i''» WjajCj = cic^cs- The sign of t 
depends on that of v, unless = 0, in which case the term vanishes. 



7.2 Kinetic orders with common influence speciflcation 

In the previous subsection we studied injectivity of a network over the set of kinetics K,g{J\f)[\'] 
for a fixed kinetic order v. This family contains all general mass- action kinetics with kinetic 
order v and arbitrary rate vector k. The results therefore assert injectivity independently of 
the rate constants. In this subsection we are concerned with asserting injectivity for sets of 
general mass-action kinetics determined by influence specifications associated to kinetic orders. 
Specifically, we will characterize networks that are injective over any set of general mass- 
action kinetics for which the kinetic orders share common associated influence specification. 
Furthermore, we will use this characterization to provide a similar result for networks with 
associated influence speciflcations that are bounded by two given influence speciflcations. 

To proceed we need some additional notation and definitions. We first introduce a partial 
order on the set of influence specifications for a network N . One infiuence specification X is 
said to be smaller than another influence specification T if ly^yi ^ ^y^y' ^^"^ ^y^y' — ^y^y' 
for all y ^ y' € TZ. If this is the case then we write X <X. The minimal element in this order 
is the zero influence speciflcation, that is the influence speciflcation that takes value for all 
reactions on all species. There is not a unique maximal element in this order but all maximal 
elements must fulflll ly^yi U I'^yi = {1, . . . , n}. 

For amxk matrix V, let be the sign pattern of V, that is, the mxk matrix indicating 
whether the entries of V are positive, negative or zero, respectively. If the entries with positive 
(resp. negative) signs of a matrix V' form a subset of the entries with positive (resp. negative) 
signs of V then we say that 'S(V') is smaller than T,{V) and write Ti(y') ^ S(y). Consider the 
set of matrices 

= {V, an mxk matrix | T.{V') = T.{V)}. 
This set is called the qualitative class as 

Definition 7.7. An m x m square matrix V has a signed- determinant if the function 
(5 := signodet: Sy ^{-1,0, 1}, V sign(det(y')) 



22 



is constant. If the image of 5 is not zero then V is called sign-nonsingular (SNS). 

In other words, any matrix with the same sign pattern as a matrix V with signed-determinant 
has a determinant of the same sign as det{V). For example a 2 x 2 matrix with one negative 
entry and three positive is SNS, whereas one with all four entries positive is not. If one column 
of a matrix is identically zero, then the matrix has signed-determinant but is not SNS. 

Any matrix V' in Ey is described by the specification of the absolute values aij = Wi j\ for 
each non-zero entry of V, such that v'^j = sign(fjj)ajj. This gives a positive parameterization 
of Sy. The function det(-) restricted to Sy becomes a polynomial function in the variables ajj 
for Vij 7^ 0, which we denote by pv{aij). In fact, if 1/ is an m x m matrix, py(aij) is either 
the zero polynomial or a homogeneous polynomial of degree m in ajj. Further, the degree of 
each monomial in each variable Uij is either zero or one. 

If V is any matrix in Sy, then det(y) = pyd^i This equality remains valid for matrices 
V such that T,{V') ^ '^{V), now with some variables Uij vanishing. 

Lemma 7.8. Let V be an m x m square matrix. Then: 

(i) V has a signed-determinant if and only if pv is either the zero polynomial or the non-zero 
coefficients of the monomials of py have sign 5{V). 

(a) Let V be an m X m matrix such that T,{V') ^ '^(V). 

(a) IfV has a signed-determinant then V' has a signed-determinant. 

(b) IfV is SNS and S(y') = then V is SNS. 

(c) IfV and V are both SNS then d{V') = 6{V). 

(d) IfV is SNS and V has a signed-determinant then V is SNS. 

Remark 7.9. If we consider a kinetic order v = {vy^yi)y^yi^Ti to be a matrix with columns 
Vy^yi then the sign pattern is the same for all kinetic orders v' with X(v') = X(v) and, further, 
E(v') ^ S(v) if and only if X(v') < 2^(v). We will be interested in the matrices Z{R,v)j for 
some sets R and J. The columns of Z{R,v)j agree with a subset of the columns of v, hence 

S(Z(i?,v)j) = {Z{Ry)j\ X(v') =X(v)}. 

Also, E(v') ^ S(v) implies that S(Z(i?,v')) ^ T.{Z{R,v)). 

Note that a network M is injective over K.g{J\f ,T) if and only if is injective over /Cg(AA)[v] 
for all kinetic orders with X(v) = X (Lemma |6.4[). Hence, we can use Proposition |7.5| and 



Lemma 7.8 to derive a determinant criterion valid for fCg{J\f,I). Similarly, if Xi,X2 are influence 
specifications for J\f such that Xi < I2 then we can use Proposition |7.5| and Lemma [7.8| to derive 
a determinant criterion for the network to be injective over Ux|Xi-<x-<X2 ^gi-^^-^)- Lemma 17. 8| 
is applied to the matrices V = Z{R, v) j and V = Z{R, w) j with kinetic orders v, w such that 
X(v) = X and X(w) ^ X(v), respectively, in the two cases. 
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Proposition 7.10. Let N = {S,C, TZ) be a network with influence specification I and choose 
a kinetic order v with X(v) = I. Then N is injective over ICg{M,I) if and only if the following 
two statements hold: 

(i) M is injective over ICg{J\f)[v]. 

(a) For all sets R of s reactions from TZ and J £ Od{M), ifT{R)j is non-singular then 
Z{R,v)j has a signed- determinant. 

In particular the proposition is true by choosing the kinetic order v = e(X). The proposition 
provides a characterization of injectivity of a network N over ICg{M,I). The first condition 
guarantees that all coefficients in the polynomial expression of the determinant det( Jc(/k,v)) 
have the same sign or are zero, and that at least one coefficient is non-zero. Using Lemma [7.8[ ii), 
the second condition ensures that this property is preserved for all kinetic orders v with X(v) = 
I. 



Remark 7.11. Note that by virtue of Theorem 7.1 we have the following. If the determinant 
det( Jc(/k,v)) does not vanish for all (k,v) G /Cq ^ ICg{M,I) then M is injective over /Co, 
irrespectively whether any of the matrices Z{R,v)j has a signed-determinant. For example, 
if det( Jc(/k,v)) is not zero provided that Vy^y'^i > then J\f is injective over /Co = 

Proposition 7.12. Let M = {S,C,TZ) be a network and Xi ^ I2 two influence specifications 
for N . Choose kinetic orders vi,V2 with = X\ andXiy^) =^2- 

Then J\f is injective over Ux|Xi-<x^X2 ^g(-^'-^) ^/ ^'^'^ '^'^^2/ ^/ the following conditions hold: 

(i) M is injective over lCg{N')[vi] and ICg{M)[v2]- 

(a) For all sets R of s reactions from TZ and J £ Odi-M), if T[R)j is non-singular then 
Z{R,V2)j has a signed-determinant. 

In particular the proposition is true for the kinetic orders vi = e(Xi) and V2 = e(X2). Propo- 
sition 7.12 provides a characterization of injectivity of a network M over Ux|Xi-<x-<X2 ^g{-^)^)- 
The second part of the first condition guarantees that all coefficients in the polynomi al ex - 



pression of the determinant det( Jc(/k,v2)) have the same sign or are zero (Proposition 7.5). 
By Lemma |7.8[ ii), the second condition ensures that this property is preserved for all kinetic 
orders w with X(w) ^ X2 . The first part of the first condition ensures that at least one term is 
non-zero for all w with Xi ^ X{'w). Thus, the first condition could be replaced by the condition 
that there is at least one non-zero term for some kinetics (k, vi) G K.g{J\f ,Xi). Consequently, 
Proposition |7.12 can be rephrased to obtain: 



Corollary 7.13. Let N = {S,C,TZ) be a network andXi,X2 two influence specifications for J\f 
such thatXi ^ X2. The following statements are equivalent: 
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(i) M is injective over \Jx\ii^i^i2 ^g(-^ ■ 

(a) There is a £ {0, 1} such that (—1)'^ det{Z{R, V2) j) det(r(i?)j) > for all sets R of s reac- 
tions from TZ, J £ Od{M) and V2 with I{v2) =I2- Further det{Z{R,vi)j) det{T{R)j) ^ 
for at least one set R, some J E Od{M) and some vi with I{vi) = X\. 

Propositions 7.10 and |7.12" provide the following corollary. 

Corollary 7.14. Let M = {S,C,TZ) be a network and Ii ^ X2 two influence specifications for 
N . The network N is injective over }Cg{J\f,Ii) and ICg{J\f, I2) if and only if N is injective over 
Ux|x,^x^x,/Cg(AA,X). 

A natural choice for the smaller influence specification Xi is in many contexts the complex 
dependent influence specification Xq given by the kinetic order y. Proposition 7.10| and the 
discussion above also imply the following corollary. 

Corollary 7.15. Let N he a network. Assume that N is not injective over the set of mass- 
action kinetics ICg{Af)[y] and that the determinant det(Jc(/K,y)) is not identically zero. Then 
N is not injective over /Cg(AA)[v] for any kinetic order v such that Xq -< X(v). 

Remark 7.16. Assume that 0{N) = {l,...,n}, that is, the network N has no outflow 
reactions. Then the determinants of Z[R, v) j and T{R)j for all reactions sets R and all subsets 
J C {1, . . . , n} of cardinality s are exactly the s x s minors of ZiJZ, v) and r(7^), respectively. 



Example 7.17. Consider the futile cycle with the kinetic order v introduced in Example 6.5 
The kinetic orders w = (wi)i=i,...,6 with X(w) = X(v) are of the form 

wi = {ai,s\gn{v)w, 02, 0, 0, 0), W2 = (0, 0, 0, 0, as, 0), lus = (0, 0, 0, 0, 04, 0), 
= (0, 05, 0, ae, 0, 0), = (0, 0, 0, 0, 0, 07), wq = (0, 0, 0, 0, 0, ag). 

with w > Q and Oj > 0, i = 1, . . . , 8. The kinetic orders w with Xq < X(w) < X(v) also include 
the possibility w = Q. The network has no outflow reactions. Hence, using Remark 7.16 



the products det(Z(i?, w) j) det(r(ii)j) are the products of the corresponding 3x3 minors of 
Z(7^, w) and ^(7^). We have 
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The products of the corresponding 3x3 minors of the two matrices are: 

(0, — a2a4a5, — sign(T;)wa4a6, — a2a4a6, — a2a4a7, — a2a4a8, —aia^ag, —a2aQa8, —asaQOg, —a^aQOs) 
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If u > then the non-zero coefficients have negative sign. In that case it follows from Proposi- 
tion 7.10 and Lemma [7.8[ i) that the futile cycle is injective over lCg{M ,T{w)) and in particular 
over ICg{J\f,Ic) (corresponding to t; = 0). If t> < then the term — sign{v)wa4aQ is positive 
while the rest are negative and the futile cycle is not injective over /Cg(AA,X(v)). 

Consider now the kinetic order w with Wi = for i 7^ 1 and wi = (1, —wn, 1, —W12, 0, 0) 
with wii > 0. That is, we modify the kinetic order v above to incorporate a negative influence of 
5*4 in the first reaction. Let I = X(w) . We repeat the procedure above and consider the family 
of kinetic orders with associated influence specification I. The products of the corresponding 
3x3 minors have the form: 



(0, —020405, — t(;i20405 -I- 1^110406, —020406, —020407, -^120407, -020403, —Wi2a4as, 

- oioeag, -a2«6a8, -asanas, -040603). 



If a kinetic order w' fulfills I 



X(w') and ^1205 > wuae, then all terms are negative and 

the 



thus is injective over /Cg(AA)[w'] (Proposition 7.5 see also Remark 7.11). However 
term — i(;i20405 -|- u;ii0406 shows that one of the matrices Z{R,w)j does not have a signed' 



determinant (Lemma 7.8 ^i)). Therefore, M is not injective over }Cg{J\f,X). 



Example 7.18. Consider the futile cycle with the reaction dependent influence specification 
Xji. Any kinetic order v such that X(v) = X-j^ provides 
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with Vi > 0. Considering r(7^) as shown in the previous example, the non-zero products of 
corresponding 3x3 minors of these matrices have all negative sign and there is at least one 



that is non- negative. It follows that Proposition 7.10 holds and that the futile cycle is injective 
over ICg{M,Xn). The previous example showed that the futile cycle is injective over K.g(J\f ,Xc)- 
Hence, it follows from Corollary 7.14 that the cycle is injective over \^x\Xc^x-<Xtz ^giN 



The Mathematica code implementing the algorithm to decide whether the futile cycle is 
injective over ]Cg{Af,XTi) is shown in Figure [3j 



Remark 7.19. In general a network will not be injective over the set of all general mass- 
action kinetics JCg{N). In the case of the futile cycle (Example |7.6[ ), the term t = vcid^c^kik^k^^ 
is the only term depending on z;. It changes sign with v whereas none of the other non-zero 
terms do. Hence, the futile cycle is not injective over K,g{M). 
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Define the stoichiometric matrix r(7^): 



R = Traiispose[{{-l, 0,-1, 0,1,0}, {1,0, 1,0, -1,0}, {1,0, 0,1, -1,0}, 

{0,-1, 0,-1, 0,1}, {0,1, 0,1, 0,-1}, {0,1, 1,0, 0,-1}}] ; 
{s,n} = {MatrixRankCR] , Length [R[[l]]]}; 

Define the matrix Z{TZ,v): 

Infl = Transpose [{{v[l] ,0,v[2] ,0,-v[3] ,0},{-v[4] ,0,-v[5] ,0,v[6] ,0}, 
{-v[7] ,0,0,-v[8] ,v[9] ,0},{0,v[10] ,0,v[ll] ,0,-v[12]}, 
{0,-v[13] ,0,-v[14] ,0,v[15]},{0,-v[16] ,-v[17] ,0,0,v[18] }}] ; 

lengthv = 18; 

Compute the products of the corresponding minors (only non- repeated ones are kept): 

{detZ,detG} = {Minors [Infl, s] , Minors [R,s]}; 
coeffs = DeleteDuplicates [Flatten [detZ*detG] ] ; 

Check the signs of the coefficients: 

Rules = {}; 

For[i=l,i<=lengthv,i++,AppendTo [Rules, v[i] ->1] ;] ; 
sign = DeleteCases [DeleteDuplicates [Sign [Flatten 

[MonomialList [coeffs]] /.Rules]] ,0] ; 
If [Length [sign] == 1, 

Print ["The network IS injective"] , 

Print ["The network is NOT injective"]; ]; 

Figure 3: Mathematica code implementing the algorithm to decide whether the futile cycle is 
injective over K,g{N iX-j^). 

7.3 Interaction networks defined by the set R of reactions 

The conditions presented in the propositions in the previous sections relate to subnetworks of 
the network M = {S,C,7l). Let R he a set of s reactions from TZ and let Mr = {Sr,Cr,TIji) 
be the network restricted to the reactions in R, that is: Sr = S, Cr = {z\ y ^ y' e R, z = 
y or z = y'} and TZr = R. We say that Mr is the restriction of A/" to R. By definition there 
might be species of Mr that are not involved in any reaction in TZr. 

Theorem 7.20. LetAf = {S,C,TZ) be a network, v a kinetic order and let Mr he the restriction 
of M to a set R of s reactions from TZ. Assume that the stoichiometric subspace of Mr has 
dimension s and that M is injective over ICg(M,v). Let vr = ivy^y')y^y'eR- Then either 

(i) Mr is injective over ICg{MR,'VR) and all steady states are non- degenerate, or 
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(ii) Mr has only degenerate steady states. 



The theorem relates to Ref. |8t Cor. 8.1, Cor. 8.2], where injectivity of a network is studied 
relatively to injectivity of the network augmented with the "missing" outflow reactions. Also, 
the theorem relates to Ref. [Uj, where it is shown that if a restricted network Mr with stoi- 
chiometric dimension s has multiple non-degenerate (positive) steady states then so does M. 



However, Theorem 7.20 cannot be used to draw the same conclusion since non-injectivity does 



not imply that there are multiple steady states. 

Example 7.21. Consider Example |4.9[ The stoichiometric space has dimension one and there 
are two possible choices of reaction sets, Ri = {S ^ 25} and R2 = {S 0}, both with empty 
index set, Ji = J2 = 0. The networks Mr-^ and Mr2 are both injective over /Cg(A/ij, v/j), but 
is not injective over ICg{M,v). 

8 Injectivity for strictly monotonic kinetics 

In this section we extend the results on injectivity for general mass-action kinetics to cover 
X-injectivity of a network M over the set of strictly monotonic kinetics IC{M,I). 

Theorem 8.1. Let M = {S,C, TZ) be a network with influence specification X . The following 

three statements are equivalent: 

(i) M is X -injective over IC{M,X). 

(ii) M is X -injective over ICd{M,X). 

(Hi) M is injective over lCg{M,X). 

The theorem implies that for a network to be X-injective over IC{M,X) it is sufficient to 
be injective over /Co with ICg{M,X) C /Cg ^ IC{M,X). In Ref. [81 Prop. 5.2] it is shown that 
injectivity and X-injectivity are equivalent notions for mass-action kinetics. Theorem |8.1| also 
implies that Ref. [8, Prop. 5.2] holds generally, namely that M is X-injective over fC{M,X) if 
and only if is injective over IC{M,X). 

Theorem 8.2. Let M = {S,C,TZ) be a network. Then the following are equivalent: 

(i) ker(Jc(/K)) n T = {0} for all c G M!^ and K G ICd{M,X). 

(ii) ker(Jc(/K,v)) n P = {0} for all c G M*^ and (k, v) G fCg{M,X). 

If either of the two statements holds then M is X-injective over }C{M,X). 



Recall from Proposition 5.4 that ker( JcC/ii")) n P is zero if and only if det(Jc(/A')) 7^ 0. 
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Corollary 8.3. If M is injective over ICd{M ,1) then M cannot have positive degenerate steady 
states. 



The corollary follows immediately from Proposition 5.4 It also follows from Ref. [23j and 
from Ref. [8] for mass-action kinetics. 



Example 8.4. According to Example 



7.18 



1 the futile cycle is injective over Uxix^^x^Xtj ^gi-^^-^)- 
It follows from Theorem 8.1 that the tutile cycle is X-injective over fC{J\f,I) for any X such 
that Xc ^ X ^ I-ji. That is, for any kinetics that is strictly monotonic with respect to an 
influence specification that is positive exactly for the species in the reactant complexes and 
may be negative for the species in the product complexes, the futile cycle is X-injective and 
cannot admit multiple positive or non-overlapping steady states. 



9 Extensions to other types of influence speciflcations 



In Ref. [23] the authors introduce the term "weakly monotonic kinetics" which in some sense 
imposes a weaker requirement on the kinetics than the term "strictly monotonic kinetics" 



introduced here (Definition 4.6 and Lemma 4.8). In this section we assume that X is an 



influence specification such that supp(y) C Iy_^y, U Iy_^yi (see also Remark 4.2). 



Definition 9.1 (Ref. p3]). A kinetics K = {Ky^yi)y^yi^n for a network N = {S,C, TZ) is 
weakly monotonic with respect to an influence specification X if, for each pair of non-overlapping 
concentration vectors a,b £ Qk, the following implications hold for all reactions y ^ y': 

(i) if Ky^yt{a) > Ky^yi{b) thcu sign(aj — 6j) = Iy^y/{Si) ^ for some i. 

(ii) if Ky_^yi{a) = Ky^yi{h) then Oj = hi for all i E supp(y), or sign(aj — bi) = Iy^yi{Si) ^ 
and sign(aj — hj) = —Iy^y'{Sj) ^ for some distinct i,j. 

Let ICuj (M, X) denote the set of kinetics that are weakly monotonic with respect to X. 



Using the characterization of strictly monotonic kinetics provided in Lemma 4.8, we find 



that the two definitions differ in (ii) , where it is required that Oi = bi for all i € Iy_^yi U ly^y' 



and not just for i G supp(y). In this sense. Definition 9.1 imposes a weaker requirement on the 
kinetics than Definition 14.61 and 

/C(AA,X) c /C^(AA,X). 

Note that in Ref. [23], an infiuence specification X must fulfill Iy^yi{Si) = 1 for i G supp(y), 
which is not required in this paper. Definition |4.6| stipulates that all species play an equal role 



in the definition, whereas Definition 9.1 singles out the species in the reactant complex to have 
special importance. Below we show that our determinant criterion also applies to the broader 
definition of infiuence specification. 
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The determinant criterion derived in Corollary |7. 13| can be adapted to derive a determinant 
criterion for a network N to be X-injective over /C^(A/',X). We first note that the influence 
specification given by Iy^yi{Si) = Iy^yi{Si) for i G supp(y) and zero otherwise is a minimal 
element among all influence specifications for M that are smaller than X. Also note that if a, 6 
are X-non-overlapping then they are X'-non-overlapping for all influence specifications X' such 
that X' < X. Hence, we have that 

X' <X implies IC^{M,X') Q ICu,iM,X). (9.2) 

Lemma 9.3. If K & /C^(A/',X) is a general mass-action kinetics then there isX <X' ^X such 
that K e )Cg{Af,X'). That is, 

IC^ {Af, X) n Kg {M) = Kg {M, X'). 

T\X<I'<E 

Theorem 9.4. Let M = {S,C,TZ) be a network with influence specification X . The following 

statements are equivalent: 

(i) M is X-injective over K,yj{M ,X) . 
(a) M is X-injective over Uj-,^j_,j-,_,j-IC{M ,X') . 
(Hi) M is injective over Uj-,^j-^j-,^j-ICg{M ,X') . 



Together with Theorem 7.13 we derive a determinant criterion for a network M to be X- 
injective over /C^(AA,X). Further, it is straightforward to derive statements similar to those in 
Section 1] for )C^{Af,X). 



Remark 9.5. In Ref. |23j the authors introduce the concepts of a concordant network and 
of a network being concordant with respect to an influence specification X. To be concordant 
depends on the kinetics associated with the network only through the influence specification 
X. Ref. [23] shows that to be concordant with respect to X is equivalent to be injective over 



KLw{M,X). Theorem 9.4 and Theorem 7.13 provide an equivalent characterization in terms of 
the influence specification through properties of the matrices Z{R,v)j and r{R)j and bring 
out an explicit relationship to the set of general mass-action kinetics. 



10 The P-matrix property 

In Ref. [2j, a criterion is introduced that relates to injectivity of a network M = {S,C,Tl). The 
authors study a very broad class of kinetics termed non-autocatalytic (NAC) kinetics. In our 
terminology a NAC kinetics fulfills that (a) no species appear both in the reactant and the 
product complex of a reaction and (b) it belongs to Ux^x ^{■^i'^)-, where X?^ is the reaction 
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dependent influence specification. By (a), we have that X?^ is given by ay^yi{Si) = sign(yj — y^), 
y ^ y' & TZ. We impose the same constraints in this section and consider kinetics such that (a) 
and (b) are fulfilled. This implies that the qualitative class of —Z{R, v) is in the qualitative 
class of r(i?) for any R OTZ and X(v) ^ In. 

A square matrix is said to be a P-matrix if all principle minors of the matrix are positive. 
If the principle minors are non-negative the matrix is said to be a Po-mathx. Let K he a NAC 
kinetics for M. The criterion in Ref. Th. 4.1, Th. 4.3] states that 

(*) any square submatrix of the stoichiometric matrix r(7^) is SNS or singular if 
and only if —JcHk) (minus the Jacobian) is a PQ-m.atiYx.. 

Here r(7^) is the matrix with columns y' — y for all y ^ y' ^ IZ. UN contains outflow reactions 
for all species S'j — )■ 0, i = 1, . . . , n, then (*) is equivalent to —Jcifx) being a P-matrix (Ref. [21 
Cor. 4.2, Th. 4.4]). Using the results in Ref. [9j, it follows that N is injective. If the criterion 
is fulfilled then the stoichiometric matrix is said to be strongly sign determined or SSD. 

Thus, the criteria in Ref. [2] ensures injectivity when all outflow reactions are in IZ. In 
particular, this implies that the stoichiometric subspace T has maximal dimension (equal to 
the number of species n in the network). If T does not have maximal dimension then Jcifx) 



does not have maximal dimension either (equation (3.8) ) and hence —Jc{fK) is not a P-matrix. 
Consequently, the criterion cannot be used to infer injectivity in this case. However, we show in 
the next theorem that the Jacobian of a modified version of fx has full rank and is a P-matrix 
when is a general mass-action kinetics. Let Jk be fx with the last s coordinates multiplied 
by —1, that is 

fK{c) = (o;^ ■ C, . . . ,Uj'^ ■ C,-fK,d+l{c), . . . ,-fK,n{(^))- 

It follows that the Jacobian Jcifx) of fx is the Jacobian JcHk) with the bottom s rows 
multiplied by —1. 

Theorem 10.1. Let M = {S,C, TZ) he a network and X an influence specification. Let s be 
the dimension of the stoichiometric subspace and d = n — s. The following two statements are 
equivalent: 

(i) M is injective over lCg{M,I). 

(a) There exist an order of the species in S and a reduced basis o/ T"*" such that the species 
Sd+i, ■ ■ ■ , Sn appear in some reactant complexes. Further, with the associated extended 
species formation rate function f^^v, the Jacobian JcifK,v) is a P-matrix for all kinetic 
orders v with X(v) = X and k G MJ. 

We can replace lCg{J\f,X) with Ux|Xi ^x-<X2 ^9 ^'^^ -^('^) — with Xi ■< X(v) ■< X2 
in (i) and (ii), respectively, and the theorem remains valid. Under the NAC hypothesis, the 
criteria stated in Propositions |7.10| and 7.12 for arbitrary networks with general mass- action 
kinetics relate to the criterion (*): 
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Theorem 10.2. LetN = {S,C,TZ) be a network and letl-ji be the reaction dependent influence 
specification. Assume that every square submatrix ofVilZ) is SNS or singular. Let R be a set 
of s reactions from TZ and J G Od{M), as in Proposition 1 7. S\ Then: 



(i) IfV{R)j is non-singular then Z{R,v) J has a signed- determinant for all withl{v) :<Iti. 

(ii) (-l)Met(Z(i?,v)j)det(r(i?)j) >0 for allv withl{v) ^X-n. 

If one term in (ii) is positive for the kinetic order eiX-ji) then N is injective over ICg{J\f ,T-ji) 
(Proposition 7.10| ). If one term is positive for the kinetic order y then A/" is in fact injective 



over L)x\Xc-<x-<Xn^9i-^ ^-^) (Proposition 7.12). In Ref. [2] positivity is guaranteed by imposing 
the network to contain all outflow reactions. However, if this is not the case then we are only 
guaranteed that the terms are non-negative and that the Jacobian of /k,v is a Po-maiiix. 

In summary, we have shown that the criterion (*) in Ref. [2] implies our determinant 
criterion. Our criterion is weaker in that we only consider submatrices of r(7^) given by sets R 
and J, whereas Ref. |2| considers all square submatrices of r(7^). The additional requirement 
that one term is non-zero is necessary (and sufficient) to guarantee that M is injective. In 
addition our approach ensures that if M is injective then there cannot be degenerate steady 
states. Importantly, our results do not require the stoichiometric subspace to have maximal 
dimension but hold for arbitrary networks. 



11 The interaction graph 

An interaction network puts explicit constraints on the modeling equations of a dynamical 
system. We have used these constraints to derive criteria for the preclusion of multistationarity. 
However, conditions for the preclusion of multistationarity have also been given for generic 
dynamical systems described by ordinary differential equations and we will here review one 
condition that closely relates to our work.^^^ 1^ This condition is also based on the Jacobian of 
the system and takes the form of a graphical condition. For this we need some preliminaries. 

Let a dynamical system c = F{c) = (Fi(c), . . . , Fn{c)) be given with c — (ci, • . . , Cn ) G O C 
M", where is a product of open intervals of M and Fi is differentiable in the interior of Q. The 
interaction graph G(c) at c is the labeled directed graph with node set {1, . . . , n} and labels in 
the set {—1, 1} such that there is an edge from node j to node i if dFi{c) /dcj ^ 0. If this is 
the case, the edge has label given by signed Fi{c)/dcj). 

A circuit in G(c) is a sequence of distinct nodes ii, . . . ,im such that there is an edge (positive 
or negative) between i^ and ik+i for all /c < m — 1 and one between im and ii. Two circuits are 
disjoint if they do not involve the same nodes. A A;-nucleus is a collection of disjoint circuits 
which involves 1 < k < n nodes. The sign of a circuit is the product of the labels of the edges 
in the circuit. The sign of a /c-nucleus is (— 1)^"*"^ where p is the number of circuits in the 
A;-nucleus with sign equal to 1, i.e. positive.'^Sl A fc-nucleus is variable if one edge in one of the 
circuits does not have constant sign in c. 
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In Ref. [16j [Th. 2] (see below), a mild regularity condition is imposed on F. To keep the 
presentation clear, the reader is referred to the original paper for its description. We refer to 
it as condition (C). The following is true: 

Theorem 11.1 (Ref. [16J). Assume that the system c = F{c) has two non- degenerate steady 
states and that F fulfills condition (C). Then one of the following statements is true: 

(i) There exists c G 17 such that G{c) has two n-nuclei of different sign. 

(a) There is a variable n-nucleus. 

The existence of non-degenerate steady states implies that the Jacobian of F is not singular 
for all c and that a n-nucleus exists for some c. If neither (i) nor (ii) above are fulfilled then 
the system cannot have multiple non-degenerate steady states. Note that in the statement 
there is no a priori restriction to pairs of positive steady states (unless consists of positive 
vectors only) . If (ii) is not fulfilled then preclusion of multiple steady states must follow from 
the failure of (i). Therefore, we assume that G(c) does not depend on c, that is G{c) = G for 
all c (and (ii) is not fulfilled by hypothesis). 



Theorem |11.1 can then be rephrased as a statement about the preclusion of multiple non- 
degenerate steady states in any dynamical system with interaction graph G. We represent the 
interaction graph G in the form of a matrix {5j,i}j,j=i,...,n such that gj^i is the label of the edge 
from node j to node i, being zero if such an edge does not exist. 



In Theorem 11.3 below we show that failure of condition (i) in Theorem 11.1 for any such a 
graph is equivalent to a certain type of interaction networks not being injective over the class 
of differentiable kinetics. For that, we start by associating a collection of interaction networks 
with any labeled directed graph. 

Definition 11.2. Let G be a labeled directed graph with node set {1, . . . , n} and labels in the 
set {-1,1}. 

• Let J = {i\ gj^i = for all j} be the set of nodes with no incoming edges. The network 
Mg = (5, C, 7^') consists of S = {Si,..., C = 5 U {0} and 71 = {0 ^ Si, Si ^ 0\ i ^ 
J}. 

• For i ^ J, let H{i) = {j\gj^i 7^ 0}. To every partition of H{i) into two disjoint sets, 
H{i) = Hi{i) U H2{i), we associate an infiuence specification Ihj^,H2 = {Iy-^y')y^y'eTl on 
Mg by 

h^s^{Sj) = gj,i if j G Hi{i), Is^^oiSj) = -gj^i if j G H2{i), 
and zero otherwise. 
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If G is the (constant) interaction graph associated to a dynamical system as above then 
the set J corresponds to the set of indices for which Fi is constant. If J is not empty then the 



Jacobian of F is necessarily singular and Theorem 11.1 cannot be applied. By construction, 



the dimension s of the stoichiometric subspace of Mg agrees with the cardinality of J. 
Theorem 11.3. Let G be a labeled directed graph with node set {1, . . . ,n} and labels in the set 



{ — 1,1}. With the definitions of Definition 11.2 we have: 



(i) For any K G 1C(i{.Mg,^Hi,H2)j species formation rate function fx has interaction 
graph G. 

Let s be the cardinality of J . The following two statements are equivalent: 
(a) Ng is iH^^Hi-i'njective over K,d{NG,FH^^H2)- 

(Hi) G has at least one s-nucleus and all s-nuclei of G have the same sign. 

It follows that if (iii) is fulfilled then the network Mg is injective over ]Cd{NG^T) for any 
influence specification X obtained from any partition of H{i), i ^ J. 

Consider the case s = n. We have shown that if Ng is X-injective for any chosen influence 



specification X, defined as in Definition 11.2, then any dynamical system c = F{c) with associ- 
ated interaction graph G cannot have multiple non-degenerate steady states. Indeed, if Mg is 
X-injective then Theorem 11.3| ^iii) is fulfilled and consequently. Theorem ll.l| ^i) does not hold 



(assuming that F fulfills condition (C)). 

If we chose the partition with Hi[i) = H{i) and H2{i) = then any kinetic order v with 
X(v) = X is identically zero on the reactions 5^ — )• 0. Consequently, all s x s minors of Z{R, v) 
vanish whenever R includes an outflow reaction Si — )• 0. Therefore, by Propositions |7.5| |7.10[ 



and Theorem 8.1, Mg is injective over ICd{MG,I) if and only if the network Mq, defined as 
the restriction of Mg to the set of inflow reactions S'j — t- 0, i ^ J, is injective over the induced 
influence specification X'. 

Therefore, we have proved the following corollary: 

Corollary 11.4. Let M be the network consisting of n inflow reactions ^ Si and let I be an 

influence specification for M . Further, let c = F(c) be any dynamical system in M"" such that F 
fulfills condition (C) and has constant interaction graph G = {gj.i} with gj^i = lo^SiiSj). If M 
is injective over ICg{M,I), then c = F{c) cannot have multiple non- degenerate steady states. 



Since the dimension of the stoichiometric space of M in Corollary 11.4 is n, there is only 
one set of n reactions. By the definition of the infiuence specification, any kinetic order with 
associated influence X has the signs in the columns of the matrix G. Therefore, M is X- 
injective over lCg{M,X) if and only if the matrix associated with the interaction graph G is 
SNS. Therefore, 
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Corollary 11.5. Let G be an n x n SNS matrix. Then, any dynamical system c = F{c) in 
that fulfills condition (C) and has constant interaction graph G = {gj^i} cannot have multiple 
non- degenerate steady states. 

Importantly, the corollary is saying that injectivity over the class of general mass-action 
kinetics suffices to preclude multistationarity in arbitrary dynamical systems, not only those 
originating from the specification of an interaction network and a kinetics. A possible extension 
to the theorem considers all influence specifications that are smaller than the one defined by G 
and, hence, it would assert injectivity with respect to /C^(A/',X). 

Example 11.6. Consider Example 4.4 In Ref. [15] the gene network depicted in Figure [2] is 
modeled as 

ai . . a2Ci . a3CiC2 ^ 

I + P1C3 1 + /32C1 (1 + P3Cl)(l + /34C2) 

for positive parameters, ai,/3i,5i. Here degradation of each gene Gi,G2,G3 is incorporated. 
The interaction graph G associated with the system is constant with matrix 




This matrix is SNS and therefore, by Corollary 11.5, the system cannot have multiple non- 
degenerate steady states. 

Alternatively, the model arises as the dynamical system associated with the interaction 
network with reactions — )• Gi,Gi — )• 0, influence specification I = [Iy^yi)y^yi^'ji with (only 
non-zero values are listed) 

ICi^oiGi) = 1 for all i , 10^01(^3) = —1, /o->G2(G'i) = lo^G^iGi) = Iq^G'AG2) = 1, 
and kinetics K = {Ky^y/)y^y' fzn, with Kq^^q = 6iCi and 

I + PIC3 I + P2CI (1 -h /33Cl)(l + /34C2) 



We can use Proposition 7.10 and Theorem 8.1 to conclude that the interaction network is 
injective. 

Note that in view of Corollary 11.4, the first procedure corresponds to analyze injectivity 
of the network with only inflow reactions and influence defined by the entries of G. Since 
a kinetics must be positive, the entry —1 in position (1, 1) of G corresponds to a decreasing 
kinetics in ci for the inflow reaction — )• Gi. This is different from the system we started from, 
which had a negative summand (— (5iCi). In other words, the dynamical systems characterized 
in Corollary [11.4 have positive species formation rate functions. 
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As we can deduce from the results above, preclusion of multistationarity by the methods 
of Ref. [in] is essentially preclusion of multistationarity in networks in which only inflow reac- 
tions are considered. Knowledge about the underlying network structure allows us to preclude 
multistationarity for a bigger class of dynamical systems. We do not only "see" the signs of the 
entries of the Jacobian, but also the terms that contribute to the signs. That is, if the network 
has a reaction of the form A + B ^ C then the species formation rate function has the term 
Ksi+S2^S3ic) in the entries 1,2,3 for any kinetics and it is negative for the first two. This is 
illustrated in the following simple example. 

Example 11.7. Consider the interaction graph G with associated matrix ^ ^"'^ ^1 ) ' '^^^^ 

matrix is not SNS and hence preclusion of multistationarity cannot be achieved by the methods 
in Ref. [I6]. Consider now the network J\f with set of reactions TZ = {Si — )• 5*2 , 52 — )• 0} . The 
dimension of the stoichiometric subspace is n = 2. For any kinetics K, the ODE system 
associated with the network is of the from 



Consider the influence specification with non-zero terms Is2^o{S2) = 1, Isi-^S2iSi) = 1 and 
l5'^_5.5'2 ('S'2) = —1- li K £ lCd{M,I) then the Jacobian of the species formation rate function 
has interaction graph G. The only set of reactions of cardinality n = 2 is TZ itself. For any 
kinetic order v such that I = X(v), the matri x Z{ TZ, v) has sign 1 s ince the sign pattern is 
1 ' 



-1 1 



Thus, it follows from Proposition 



K(i{M,I) and multistationarity cannot occur. 



7.10 



and Theorem 



.1 



that M is injective over 



12 Hill-type kinetics and injectivity 

Let M = {S,C, TZ) be a network and I an influence specification for A^. The key to the 
statements in Section |8] is that whenever there are two non-overlapping concentration vectors 
a, b such that a — 6 G F and fxio) = fxib) for some kinetics K G /C(AA, X), then we can find two 
positive concentration vectors d,b such that a — b £ T and fK,v{a) = /k,v(&) for some general 
mass-action kinetics {k, v) G lCg{M,I). However, the latter property could be fulfilled by many 
classes of kinetics other than the class of general mass-action kinetics. 

One such class of kinetics is Hill-type kinetics. We say that a kinetics is of Hill-type with 
respect to an influence specification Z ii K = (Ky^yi)y^yi^ii takes the form 

Ky^y'{c) = ky^yi JJ - ^y^y',v ' 

i=i "y^y',i + 
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with c G M_|_ (defined by continuity at the boundary), ky^yi G M+, 5y^yi G M_|_ and Vy^yi G M" 
for all y ^ y' £ TZ, such that 

SUPP^ {Vy^y>) = ly^y,, SUpp'{Vy^y,) = ly^y,, and S\Xpp{6y^y,) = S\ipp{v y^y,) . 

The definition is very similar to that of general mass-action kinetics with the only difference 
being the factor 5y^y'^i in the denominators. A term with Vy^yi^i > defines a positive influence 
on the reaction y ^ y' while a term with Vy^y'^i < defines a negative infiuence. Compared 
to general mass-action kinetics the constant Sy^yi^i moderates a negative influence of a species 
Si when it is in low concentration. 

Let AC = {ky^y>)y^y>(ZTl, Cl = i6y^y')y^y'e'R., SUd V = {Vy^y')y^y>(ZTl. WB dCUOte & ffill-type 

kinetics by = (k, d, v), the set of Hill-type kinetics for J\f with respect to X by /Ch(AA,X) and 
the set of all Hill- type kinetics for M by ICh{M)- Hill- type kinetics include Michaelis-Menten 
kinetics as a special case when Vy^yi^i is an integer.'^ In contrast, general mass-action kinetics 
are not of Hill-type. However, general mass-action kinetics can be obtained as a limiting case 
of Hill-type kinetics by letting ky^yi and the non-zero entries of Vy^yi tend to inflnity such 
that ky^y' /Yl^Vy-^y'^i couvcrgcs to a positive constant. 

Hill-type kinetics might be considered biochemically more reasonable than general mass- 
action kinetics as they are deflned for all M, in contrast to general mass-action kinetics that 
might not be defined for points on the boundary of M^. In addition. Hill- type kinetics or 
Michaelis-Menten kinetics are often obtained when variables (species) are eliminated from the 
modelling equations.'^ 

Theorem 12.1. Let M = (5,C,7^) be a network and a,6 G M+. Then: 

(i) For every Hill-type kinetics K = (k, d,v) G ICh{M) there exists a general mass-action 
kinetics (A,w) G K,g{J\f) such that S(v) = S(w), fK{o) = f\,w{o) and fxib) = /a,w(^)- 

(a) For every general mass-action kinetics (A,w) G ICg{M) there exists a Hill-type kinetics 
K = (K,d, v) G Kni-M) such that S(v) = i;(w), /x(a) = fx,v^{a) and fxib) = /a,w(^)- 

In particular, /A,w(a) = /a,w(^) for som,e general mass-action kinetics (A,w) G ICg{M) if 
and only if //<(a) = fxib) for some Hill-type kinetics K = (k, d, v) G K.h{-I^) such that 

i;(v) = s(w). 

As a consequence, injectivity of a network over /C(A^,X) is guaranteed by injectivity of 
M over /Ci^(A/',X). It also follows that M is injective over K.g{J\f ,X) if and only if is injective 
over /Cf/(A/',X). Furthermore, we have that N has multiple positive steady states in some 
stoichiometric class with respect to a Hill-type kinetics if and only if N has multiple positive 
steady states in the same stoichiometric class with respect to a general mass-action kinetics. 
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A Proofs 



Proof of Lemma |4.8| , We first prove the forward implication, that is, that a strictly monotonic 
kinetics fulfills (i) and (ii). For two vectors u,v in and i = 0, . . . ,n, let u (Bi v denote the 
vector consisting of the first i coordinates of u and the last n—i coordinates of v, with u(Bov = v 
and u (Bn v = u. Assume that the kinetics is strictly monotonic with respect to I. Let a, b be 
two non-overlapping vectors and consider Ky^yi{a) and Ky^yi{h) for some reaction y — )• y' . 
We have 



n 



Ky^yl{a) = Ky^yl{h) + } ^ K y^y> {t <B j- 1 «) " Ky^yl{b®i O) . (A.l) 

i=l 

If Ky^yi{a) > Ky^yi{h) thcu Ky^yi{h®i-i «) " Ky^yi{b ®i > for soffic i and in particular 
Ky^yi{b ©i-1 a) > 0. Since K respects X, it follows that Iy_^yr ^ supp(5 (Bi-i a) and hence 
6©j_ia belongs to Qxily^y/)- The two vectors 6©i_ia and b(Bia differ only in the i-th entries, 
which are Cj, 6j, respectively. If also 6©i a belongs to ^Kily_^yi), then by Definition 4.6 'iii) we 

have i S I^.„i U If i G then Ky^y/ is increasing in the i-th coordinate and hence 

y y ^y y ^y ^ ^ 

1 = sign(aj — bi) = Iy^y'{Si). It i £ ^y^y'i then Ky^yi is decreasing in the i-th. coordinate and 
— 1 = sign(aj — bi) = Iy^yi{Si). Hence (i) holds in this case. If 6©j a is not in ^Ixily^y/) then 
necessarily 6j = 0, aj > and i e ^y-^y'- Hence 1 = sign(aj) = Iy^y'{Si) and (i) holds as well. 

Assume that Ky^yi{a) = Ky^yi {b). Since a, 5 are non-overlapping, if Oj = for some 
i € ly^yn then bi ^ Q for all i G ^y^y'- Since K respects the infiuence specification X, we 
obtain that Ky^yi{a) = and Ky^y/{b) ^ contradicting the equality. Therefore, Oj / and 
bi^O for aU i G ly^yi, and hence 1^^^, C supp(a), Jjl^„, C supp(6), and C supp(6©j a) 

for all i. Consider the sum of differences in equation ( A.l ). If there is at least one positive and 
one negative term then we proceed as for Ky^yi{a) > Ky^yi{b) and conclude that (ii) holds. 
Now assume that Ky^y'{b (Bi-i a) = Ky ^yi [b (Bi a) for all i = 1, . . . , n. In particular, taking 

i G , , U , , it follows from Definition 

y^y y^y 

part of the proof. 

To prove the reverse implication, assume that (i) and (ii) are fulfilled. Let c, (i G ^K{Iy^yi) 
be two vectors that differ only in the z-th coordinate. If Ky_^yi{c) > Ky..^y'{d), then by (i) we 
have sign(ci — di) = Iy^yi{Si) ^ 0. It follows that if i G ^y^y'^ that is, Iy^yi{Si) = 0, then 
Ky^yi{-) is constant in the i-th coordinate. If i G ^y^y' ^y^y' ^'^^ Ky^yi{c) = Ky^yi{d) 
then according to (ii) we have q = di contradicting q ^ di (the second option cannot occur 
since c and d differ in exactly one coordinate). Therefore Ky^yi{c) ^ Ky^y/{d). Using (i) we 
conclude that Ky^yi{-) is increasing/decreasing in the i-th coordinate depending on the sign 
of Iy^yi{Si). This completes the proof. □ 



4.6 



that Oj = bi and (ii) holds. It completes the first 



Proof of Theorem 7.1 The equivalence of (ii) and (iii) is in Proposition 5.4 The equivalence 
of (i) and (ii) can be proved similarly to Theorem 5.6 in Ref. [8j: M is not injective over 
ICg{M)[v] if and only if there exists a rate vector k and distinct a,b £ M" such that a — b gT 
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and fn^via) = fn^vib). Further, ker( Jc(/,,_v)) n T / {0} for some c G M" and a rate vector rj 
if and only if there exists 7 G T such that Jc{fri,v)i^) = 0. By definition and using equation 



(6.2), 



Jc{fv,v){l) = <^ J2 Vy^y'c"y^'''{Vy^y' *cl){y' - y) = 0. 
y^y'eH 

It suffices to show that given distinct a,b €z M" such that 7 := a — 6 £ F and some rate vector 
K, there exist c £ and rj G such that 

ky^y' (a^'y^y' - b'"v~*v' ) = rjy^y' c'^y^y' (vy^y/ *c 7) (A. 2) 

and vice versa. Consider distinct a,b £ M" such that 7:=a — 6gF,77^0. If ai — bi ^ 
define Cj := iog'~/b ^ ^ ^^"^ Q = 1 otherwise. Since 7, = if Oj = we have that a'"y^y' = 
b'"v^y' e^y^y'*""^ . Note that the signs of a"y^y' — b"y^y' = b'"v^v' {e'"y^y'*'''^ — 1) and Vy^y' *c 7 

agree. If Vy^yi *c 7 = 0, let rjy^yi = 1. Otherwise, we let rjy^yi = ^^v^^^, -—^ — > and 

\k.2\ is fulfilled. Reciprocally, given 7 G F, 7 7^ 0, c € and r] G M?, define a, 6 G by 



6, = = 1 if 7, = and bi = ^ijie^^l'"' - \) > 0, = 6je^*/''* > otherwise. Then a-b = 7 G F 
and a^y^y' — b^y-*y' = }fy-*y' (^e'"y^y'*''^ — 1) has the same sign as y *c 7- If a"y^y' — b^y^y' = 0, 



define = 1. Otherwise, define fc,,_>,,' = — — ^7*^ ,v rr > and equality (A. 2) is 

y^y ' y^y c y^y (a y^y —b y^y ) 1 f 



fulfilled. □ 



Proof of Proposition |7.2 , The Jacobian of /k,v at c is a matrix with entries being integers 
for the top d rows and of the form ky^yic"^ for some vector m for the bottom s rows. It follows 
that det( Jc(/k,v)) is either zero or a homogeneous polynomial of degree s in the rate constants 
'^y^y ■ 

Consider now an arbitrary network N' = {S' ,C' ,TZ') with a general mass-action kinetics 
(A,w) and A = {ly^yi)y^yi^Tii . The Jacobian of /a,w takes the form 

Jc{f\^){l) = Y ly^y'c'"y^y'{wy^yi *ci){y - y) 

y^y'en 

for 7 G M". If {ei, . . . , e„} denotes the canonical basis of M", then the determinant of Jc(/a,w) 
is 

det ^ ly^yi c'"y^y' {wy^y, *c ei ) (y' - y ) , . . . , ^ ly^y' ^"y^y' {wy^yi *c en){y' - y) 

Xy^y'&TZ y^y'&l 
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By definition of the product u *c v, it follows that det {Jcifx,w)) is (recall that 1 = (1, . . . , 1)) 

c"^ det ^ ly^y>c^y^y' {wy^y> ■ei){y' -y),..., ^ ly^y,c'"v^y' {wy^yi ■ en){y' - y) , 

(A.3) 

where u ■ v is the standard scalar product in M". Each term in the determinant expansion of 



(A.3) is a scalar times a product of the form 

1=1 1=1 

for some set of n reactions R = {yi — )• y'^, . . . ,?/„ — )• y^} C 7^. Consider a fixed set R of n 
reactions and let &n denote the set of all permutations of {1, . . . , n}. Then the coefficient of 

C-^U^,ly^^y>C^y^-y^ is 

Yl ((^s/.cD^y^d) ■ - Vail)), • ■ ■ , (u'j^^j^j^j;;^^^ • e„)(2/;(„) - 

o-e6„ 



o-ee 



= Z n^J/<.{.)^!^;(0'^ (y'l - yi, . . . - yn) , 

\cre6„ \i=l / / 

det(u;j^j^j^/ , . . . , Wy^^y'J det {y[ - yi, . . . , y^ - y„) = det(Z(i?, w)) det(r(i?)), 

ni=i ^yi^y- is 



where it;,, , , ,• is the i-th coordinate of w,, , , . Thus, the coefficient of the monomial 



c ^^^^=^'"'''-"Uet(Z(i?,w))det(r(i?)), (A.4) 

Note that the argument used here is similar to the proof of Theorem 3.2 in Ref. Observe 
that if R contains repeated reactions, then det(r(i?)) is zero. Therefore det {Jc{fx,w)) is either 
zero or a homogeneous polynomial of degree n and linear in each ly^yi- 

Let us go back to the network M = (5,C,7^) with general mass-action kinetics (k,v). 
Define a network N' = {S' ,C' ,TZ') by augmenting M with the reactions in 0{M) such that 
S' = S,C' = CU {0} U {Si\ i G 0{M)] and 7^' = 7^ U {S'i ^ 0| i e 0{J\f)]. Let Wy^y' = Vy^yi 
and ly^yi = ky^yi if y — y' G 7^, and wsi^o = Cj and Isi^o = 1 if i G 0{M). Then, if 
A = {ly^y')y^y',ziii , thc palr (A, w) defines a general mass-action kinetics for A/"'. 
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With this choice of A we can view det (Jc(/a.w)) for N' as a polynomial in ky^yi, linear in 
each of the rate constant (but which is no longer homogeneous because some rate constants 
are fixed to 1). It is shown in the Proof of Cor. 7.2 in Ref. [8j for mass-action kinetics that the 
non-zero terms in the polynomial expansion of det( Jc(/re,v)) precisely are the non-zero terms 
in the polynomial expansion of det(Jc(/A,w)) of degree s and that 

(-l)'^det(Jc(/A,w)) = det( Jc(7«;,v)) + terms of total de gree smaller than s. 

The proof does not make use of the specific form of the power-law c"y^y' and thus applies 
in this context as well. It follows that det( Jc(/k,v)) is a homogeneous polynomial in the rate 
constants ky^yi of degree s and linear in each of them. 

By definition of A/"', the terms of degree s in det( Jc(/a,w)) correspond to the sets R oi n 
reactions in N' such that contain exactly d outflow reactions Sj — t- for some j G 0{N). Let 
R = {vi ~^ Hi-, - ■ ■ lUs — ^ y's} be the set of s reactions of A/" in i? and J € OdiM) be the set 
of the indices for the d outflow reactions. We can assume that the outflow reactions are the 
(ordered) reactions s + 1, . . . , n in i?. Then the corresponding summand in det( Jc(/k,v)) is 



s 

{-if n h.^yf^''^^^' "''^''^ det(Z(i?, w)) det(r(i?)) 



i=l 

s 

= ^K^y'C~^^^^='''-^^< det{Z{R,w)j)dei{T{R)j)\{c,, 

where we expanded det(Z(i?, w)) and det(r(i?)) along the j-th columns for j G J: each of the 
d reactions Sj — t- corresponds to a column in T[R) with zeros everywhere apart from the j-th 
entry which is —1. Similarly for Z{R,'w) with the only difference that the entry is 1. 

The coefficient of 11^=1 ^yi^y' the sum of all coefficients obtained from sets R containing 
R and d outfiow reactions. Therefore, the coefficient of 0^=1 ^yi~^y' det( Jc(/k,v)) is 

^-i+EU-,,^,^ ^ det(Z(i2,v)j)det(r(i?)j)nci. 

This completes the proof. □ 

Proof of Proposition 7.5[ , If (ii) holds then (i) is a consequence of Theorem 7.1 and Propo- 
sition 7.2 To show that (i) implies (ii) we do the following. We first show that for fixed c G 
the sum A{R, c) — X^jgc'^(_A/) det{Z(R, v)j) det(r(i?)j) OjeJ '^J same sign for all sets R 

of s reactions from TZ. If this is not the case we could choose rate vectors ki and H2 such that 
the sign of the determinants of Jc(/ki,v) and Jc(/k2,v) would differ. This can be done because 
the determinant of Jc(/k,v) is linear in each rate constant ky^yi (Proposition 7.2). Hence, by 
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continuity of det(Z(i?, v) j), there would also be kq such that the determinant of Jc(/ko,v) is 



zero, contradicting (i) and Theorem 7.1 



Secondly, if det(Z(i?, v) j.) det(r(i?) jj differs in sign for two sets Ji, J2 S Od{M) and 
fixed R, then we could choose concentration vectors ai and 02 in M" such that A{R,ai) and 
A{R, 02), respectively, differ in sign. This is possible because A{R, c) is linear in each coordinate 
Cj, i = 1, . . . ,n. Hence, by continuity of A{R,c), there is also b G M" such that A{R,b) = 0. 
Since A{R, b) = and (i) holds by assumption then there is R such that A{R, b) 7^ 0, say 
A{R,b) > 0. It follows that there is an open set Ub ^ around b such that A{R,-) takes 
positive and negative values and A{R, •) > 0. Taking c € C/f,, we reach a contradiction with the 
first part of the proof. □ 



Proof of Lemma 7.8, (i) V has a signed-determinant if and only if the polynomial Pvio-ij) 
has a constant sign when evaluated in positive values of aij. The equivalence follows from the 
fact that each variable has degree zero or one in pv{aij)- 

(ii) The determinant of V is given by pyd^i If V has a signed-determinant and pv = 
then V is not SNS and sign(det(y")) = for ah V" such that ^ If pv is not 

identically zero, then by (i) all non-zero monomials of pv have the same sign. It follows that 
either pyd^^ j|) vanishes or has sign S{V). Any matrix V" in the qualitative class of V' would 
fall in the same case. Statements (a)-(c) follow from this observation. If p^lbijl) does not 
vanish, then Pvi\vi,j\) 7^ as well and (d) follows. □ 

Proof of Proposition 7.10[ We will apply Lemma 7.8 to the matrices Z{R, v) j. Recall that 



X(v') = X(v) if and only if S(v') = S(v). Assume that M is injective over ICg{M,I). Then (i) is 
trivially fulfilled. To see (ii), assume that det(r(i?)j) / for some R and J. Fix V = Z{R, v) j 



and let pv be defined as in the discussion above Lemma 7.8, Let Sv = {w| E(w) = S(v)} such 
that the qualitative class of V is exactly the set of matrices Z{R, w) j for w G Ev (Remark 7.9 ). 
Let wj denote w with the entries j £ J for each Wy^yi removed. Then, det(Z(i2, w) j) = 
Pv(|wj|). We will show that pv(|wj|) has a constant sign for all w € Ev. If Pv(|wj|) does 
not have constant sign then py has at least a positive and a negative monomial. Since all 
degrees of the variables of pv are zero or one, we can by continuity choose vq € Ev such 
that pvdvo^jl) = (the set Ev inherits a Euclidean topology from the Euclidean space it is 
embedded in). Since M is injective over ICg{M,I), it is injective over ICg{M)[vo] and there must 
be some other sets Ro,Jo with det(r(i2o)jo) det(Z(i?, vq) Jq) 7^ 0, say det(Z(i?Oi vq) jg) > 
(Proposition |7.5[ ii)). In an open set fl around vq E Ev, det{Z{R, takes both positive 
and negative values. However Q can be chosen such that det(Z(i?0) v) j^) > for all v in 
Hence det(Z(i?0) v) j^) det(r(i?o) Jo) has a constant sign in J7, while det(Z(i?, v) j) det(r(i?) j) 
changes sign. We can therefore choose v G for which J\f is not injective over lCg{J\f)[v] 
(Proposition |7.5[ ii)) and we reach a contradiction. Hence (ii) is true. 

To prove the reverse we do the following. Let w be any kinetic order such that X(w) = I. By 
(ii), whenever det(r(i?)j) 7^ the matrix Z{R,v)j has a signed-determinant. Since Z{R,w)j 
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belongs to the qualitative class of Z{R,v)j, we have that 



sign{det{Z{R,w)j) det T{R)j)) = sign(det(Z(i?, v) j) det r(i?) j)). 



(A.5) 



By (i), M is injective over lCg{M)[v]. Thus by Proposition 7.5 li) all non-zero products 
det(Z(i?, v)j) detr(i?)j) have the same sign and at least one does not vanish. By ( A.5| ) the 
same holds for w and thus it follows from Proposition 7.5 that M is injective over /Cg(A/')[w] 
as well. □ 



Proof of Proposition 



7.12 



jl If is injective over ljx\Xi-<x-<i2^9(-^-'-^)-' then (i) is trivially 
fulfilled and (ii) follows irom Proposition 7.10[ ii). 



□ 



Assume that (i) and (ii) hold. If v is a kinetic order satisfying S(v) ^ S(v2) then (ii) 
together with Remark |7.9| and Lemma |7.8[ ii) imply that Z(i?, v)j has a signed-determinant 
for all reaction sets R and J G OdiM) such that det(r(ii)j) ^ 0. In particular, this state- 
ment holds for V = vi. Further, assume that v satisfies S(vi) < 5](v) < S(v2). Since 
M is injective over /Cg(A/')[vi], there exists at least one set of reactions R and J G C'd(AA) 
such that det(Z(i?, vi) j) det(r(i?)j) ^ 0. Since Z(i?,vi)j has a signed-determinant, it must 
be SNS. From Lemma |7.8| ^ii,d) it follows that Z{R,y)j is also SNS and in particular that 
det(Z(i?, v) j) det(r(i?) j) / 0. Consequently, it follows from Proposition 7.10 that M is in- 
jective over /Cg(Ar,X(v)). Since Ux|Xi^x^X2 'Cg(Ar,X) = Uv| s(vi)dS{v)dS(v2) ''^'/(■^'^(v))' 
proposition is proved. □ 

Proof of Corollary 7.15, Since the determinant det( Jc(/K,y)) is not identically zero there 
is a term in its expansion in k with positive coefficient and a term with negative coefficient. 
Since, Ic ^ ^(^)) the terms in the polynomial expansion of det( Jc(/k,v)) in k cannot have all 

implies that ICg{M)[v] is not injective. 

Assume that M is injective over /Cg(A/')[v]. Proposition 
implies that the non-zero products det(Z(i?', v) j) det(r(i?')j) have the same sign 5 for all sets 
R' of s reactions from TZ and J G Od{M), and that at least one of the products is non-zero. 

For the network Mr there is only one set of s reactions from R, namely R itself (s is 
the dimension of the stoichiometric subspace of Mr). Observe that Z{R,\-r) = Z(R,v) and 
that 0{N) C 0{Mr), but that there might be outflow reactions in TZ that are not in R. For 
J G Od{M), we have that det(Z(_R, v/j) j) d(ii{T{R)j) is either zero or has sign 5. If J ^ Odi^f) 
then there exists an outflow reaction — >• that does not belong to R but belongs to IZ. In 
this case det(r(i2)j) = it det(r(i2) j^jj}), where R is R augmented with the reaction Si — t- 0. 
Since the dimension of the stoichiometric subspace of is s and Si ^ belongs to TZ then 
—Si is a linear combination of the reaction vectors in R. Hence det(r(i?)j\{j}) = 0. It follows 
that det(Z(i2, v/j) j) det(r(i?) j) is either zero or has sign 6 for all J G Od{MR). 

If det(Z(i?, v/j)j det(r(i?)j) = for all J G Od{MR), then all steady states (if there are any) 
are degenerate (Proposition 5.4). If for some J G Od{MR) we have det(Z(i?, v/j) j det(r(i?)j) ^ 

that Mr is injective over /Cg(A/'i?)[v/j] and hence all steady 

□ 



the same sign. Thus Proposition 
Proof of Theorem 



7.20 



7.5 



0, then it follows from Proposition 
states are non-degenerate. 



7.5 
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Proof of Theorem 8.1. The proof is similar in spirit to arguments presented in Ref. |23j . 



By Lemma 4.12, (i) implies (ii), and since general mass-action kinetics are differentiable with 
respect to the influence specification, (ii) implies (iii). Let us prove that (iii) implies (i). 
Assume that M is injective over fCg(J\f,I) but not X-injective over /C(A/',X). Then there exists 
K £ K,{J\f ,X) and distinct non-overlapping vectors a,b £ M._^_ such that 7 := a — 6 G T and 



Ik {a) = fxib). Let a 



Ky^yi (a) — Ky^yi (5) for 1/ ^ i/' £ TZ . Then, 



ay^y'iy' -y)= Y {Ky^y'ia) - Ky^yi{h)){y - y) = 0. 



(A.6) 



y^y'en 



y-^y'en 



If dy^yi > for some reaction y — )• y' E 7^, then by Lemma 4.8 we have that sign(aj — h 



Iy^y'{Si) 7^ for some species Si. Similarly, if oiy^yi < then there is some species Si such 

a,- — bi =0 for all 



y-*y' 

hi) = -Iy^yt{Si) / 0. Finally, if ay^yi = then 7, 
, or there exist distinct indices i,j such that sign(ai — h. 



that sign(aj 
^ ^ ^y^y' ^ h^y 

sign(aj — hj) = —Iy^yi{Sj) ^ (here we use that o, h are non-overlapping.) 

The above consideration implies that for each y 
such that 



'y^yi{Si) / and 
y' G TZ we can find vectors Vy^y/ £ 



y^y 



supp (n 



y-^y' ) 



r 



y^y 



and a 



For example, assume that cty^y' > such that sign(a. 



y^y' 



•^y^y 



^y^y 



■ 7- 



(A.7) 



{Si) / for some i. 



If Oj > bi, choose Vy^y/j for i ^ j such that sign{vy^y'j) = Iy^y'{Sj) and /3 := Oiy^y' — 
^j^iVy^y'jiaj - hj) > 0. By defining Vy^y^^i = (3/{ai - bi) > 0, we obtain ay^y> = Vy^y, ■ 7 
and the above relations are satisfied. Similar arguments apply in the other situations. 

For all real x we have that sign(2;) = sign(e^ ~ !)• If ^y^y' "7 = 0, let r]y^yi = 1, and 
otherwise let rjy^yi = Oy^y' /{e^y^y'"^ — 1) such that 



a. 



y^y 



Jy^y' 



7 = r]y^y'{e''y^y''"' - 1) 



for all y ^ y' £ TZ. Note that i]y^y' > 0. Next we define two positive vectors a,b £ M". For 
i = 1, . . . ,n, if 7i 7^ 0, let bi be the number fulfilling 7^ = bi(e'^' ~ !)• If 7* = 0, let bi = 1. 
Define d = be"', then a — 6 = 7 G F and d/b = e"* . Note that both d and b are positive and 
from the definitions we obtain 



a. 



y-^y 



b^y^y' 



Define k = {ky^y>)y^yi( zTi b y ky^y> = r]y^yi /h^v^v' and let v = {vy^y')y^yi(zTi. Then 
(k, v) £ Kg{M,T) and from 



/K,v(a) 



^ ky^yid^y-^y' {y - y) 
y-^y'elZ 



J2 ky^y'V'y^^'iy' -y) = fnAb)- 

y^y'GlZ 
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Consequently M is not injective over ICg{M,I). 



□ 



Proof of Theorem 8.2, Clearly (i) implies (ii). Assume that (ii) holds and ker( Jc(/a')) HF ^ 
{0} for some c G M" and K G /Cd(7V,X). Then there exists c* G and 7 G T such that 
Jc* {fK){l) = 0. By equation ( |3.8[ ) and letting Vy^yi be the partial derivative of Ky^yi evaluated 
at c*, we have 



J.-Uk)(~i) 



E 



%^(c*)-7) {y'-y) 



y^y'en 



y) = 0. 



If we let K 
have 



2/' G 71 and c = (1, . . . , 1) G 

y^y' *ci){y' -y) = 0. 



we 



{ky^y')y^y>^TZ Wlth ky^y> = 1 for ah 1/ 

Jc*{fK){l)= ^ ky^y^c^'y-^y' {u 

y^y'&l 

Since the kinetics K belongs to K,d{M,X), Definition 4.11 implies that supp"^(fy^y') = ly^yi 
and supp~ (f y_s>y' ) = ly^yi- It follows that v = {vy^yi}y-^yi^Tz, is a kinetic order with X(v) = X 
and (k, v) G K.g{J\f,X). Further Jc{fK,v)il) = Jc*{fK){l) = contradicting (ii). □ 

Proof of Lemma 9.3| , Assume that K G /C^(A/', X) is a general mass-action kinetics with 
kinetic order v. Then K G }Cg{N ,X') with X' = X(v). Let us prove that X' < X. Let i be such 



that/;^^,(5i 



1, that is, Vy^yi^i > 0. Consider a, 6 G such that a,- 



>j,i j, and > 6^. 

i) 



9.1 



Then Ky^yi{a) > Ky^yi{h) because K is a. general mass-action kinetics. By Definition 
and using that a, b only diff'er in the index i, we have Iy^yi{Si) = sign(aj — hi) = 1. Therefore, 
Iy^y'{Si) = Iy^y'{Si). Wb procccd Similarly if Iy^y,{Si) = —1 to conclude that X' ^X. □ 

Proof of Theorem 9.4[ By Lemma 9.3, (i) implies (iii). (ii) and (iii) are equivalent due to 
Theorem 8.1 That (iii) implies (i) can be proved similarly to the proof of Theorem 8.1 The 
only difference lies in choosing Vy^y' G 
choose V, 



for each y ^ y' £ TZ, see equation (A.7). Now we 



"y^y 



G M" such that supp(y) C supp(?;; 



y^y ): 



and 



supp 



) c/^ 



supp {Vy-^y') C 



y^y 



and 



a. 



y^y 



'Jy^y' ■ 7- 



This can be done using an argument similar to the argument provided below equation (A.7). □ 



Pro of o f Theorem 10.1 The proof follows the idea given in Ref. [8] [Sect. 10]. By Proposi- 

M is injective over /Cg(A/')[v] if and only if det( Jc(/k,v)) for all 



tion 

cG ] 



7.5 



and Theorem 



7.1 



■+ 



, which happens if and only if 



(**) there exists a G {0, 1} such that (-1)^ det(Z(iZ, v) j) det(r(i?)j) > for all 
sets R oi s reactions from TZ and J G OdiM) and at least one of these terms is not 
zero. 
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Assume that (**) holds. Let Si-^, . . . , Si^ be the species of A/" that are not in reactant complexes 
(a complex y in a reaction y — t- y'). Then the lower s rows of Jc(/k,v) are zero in entries 

ill ■ ■ ■ lim- If (**) holds then det(o^c(/K,v)) 7^ and the columns ii,...,im of Jc(/k,v) are 
linearly independent. Since the lower s x m submatrix is identically zero, it follows that the 
upper d X m matrix has rank m (and in particular m < d). As a consequence, we can reorder 
the species in S such that Si^^, . . . , Si^ are the first m species and we are guaranteed that there 
exists a reduced basis of T-*- with that order. Note that this order is independent of the kinetic 
order v. 

Consider now the order given above. Because m < d, the species Sd+i-, ■ ■ ■ ■, Sn oi M appear 
in some reactant complex and the NAC assumption guarantees that the diagonal entries of 
Jc{fK,v) in the rows d + 1, . . . ,n are non-zero and negative. Since the diagonal entries in rows 
1, . . . , d are all equal to one, it follows that the product of the diagonal entries is a term in 
the expansion of det( Jc(/k,v)) and has sign 1 (by construction of /). The matrix Jc(/k,v) 
has full rank and is a P-matrix: if M fulfills (**) then all non-zero terms in the determinant 
expansion of Jc{fK,v) have sign 1. Each principle minor can be obtained as a sum of terms in 
the determinant expansion of Jc(/k,v) divided by the (positive) diagonal entries of the rows 
not taken in the minor. Since this sum contains at least the non-zero diagonal product term, 
all principle minors are positive. Consequently, Jc(/re,v) is a P-matrix. 

Therefore, if A/" is injective over ICg^M,!), then (**) holds for all kinetic orders v such that 
X(v) = I and consequently there exists an order of the species of M such that the matrices 
Jc{fK,v) are P-matrices. Reciprocally, if Jc{fK,v) is a P-matrix for all v then, using the results 
in Ref. 0, we conclude that /k,v is an injective function for all v and hence that M is injective 
over /Cg(AA,X). ' □ 

Proof of Theorem 10. 2[ Consider a set of reactions R, an index set J and the kinetic 
order e = e{lTi) given by Cy^yi^i = sign(yj — y^). Since every square submatrix of r(7^) is 
SNS or singular then it is in particular true for r(P)j for any R and J. The matrix Z{R,e)j 
belongs to the same qualitative class as — r(P)j. Thus, if T{R)j is non-singular, it is SNS 



and so is — r(P)j. It follows from Lemma |7.8| ^ii) and Remark 7.9 that Z{TZ,w)j has a signed- 
determinant for all w such that X(w) ^ I-ji. Further, the sign of det(Z(P, w) j) either agrees 
with the sign of det(— r(P)j) = det{Z{R,e)j) or is zero. Because r(P)j is an s x s matrix, 
we have (-l)Met(r(P) j) = det(-r(P)j), and we obtain (-l)Met(Z(P, w) j) det(r(P) j) > 
for all w such that X(w) ^ I-ji. □ 



Proof of Theorem 11.3, (i) For shortness writeX = Ihi^H2 and A/" = Mg- Let K G }Cd{N,X). 
The species formation rate function fx has entries fK,i = if i G J and fK,i{c) = i^'o^Sj(c) — 
Ksi^Q{c) if i ^ J. If i G J, then dfK,i{c)/dcj = = gj^i as desired. Assume that i ^ J . 

Since K is differentiable with respect to X, we have that sign( ^^°^^' ^ ) = I^^StiSj) and 
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sign(^^^%^^) = Is^^oiSj). By definition of X, 



dcj dcj dcj 



j i H{i) 



In each case, tlie sign is Qj^i. It follows that the species formation rate function has interaction 
graph G. 

We proceed to prove the equivalence between (ii) and (iii). The network M is X-injective 



over K,^{M if and only if M is injective over K,g{M,X) (Theorem 8.1). Hence we prove the 
equivalence between (ii) for /Cg(A/',X) and (iii). 



Recall the criteria for injectivity in Proposition 7.10 Let the dimension of the stoichiometric 
subspace of be s and note that there are s reactions Si ^ 0, i ^ J. Consequently, OdiM) = 
{J}. Consider the interaction graph G. Nodes in the set J have no incoming edges, hence 
a node i £ J cannot be part of any circuit of G. Consequently, any s-nucleus of G contains 
precisely the s nodes in {1, . . . ,n} \ J. Furthermore, an s-nucleus of G is an s-nucleus of the 
matrix G with the rows and columns with index in J removed, and vice versa. Denote this 
s X s matrix by Gj. 

Consider first the case where H2{i) = for all i. Let v be a kinetic order such that X(v) = X 
and consider Z{R,v)j for the set i? = {0 — )• Si\i J} (ordered according to i). We have that 
T(R)j is non-singular and, by construction, the sign pattern of Z(R, v) j equals the sign pattern 
of Gj. Consequently, an s-nucleus of Gj is an s-nucleus of Z{R,\-)j and vice versa. An s- 
nucleus corresponds to a non-zero term in the expansion of det{Z{R,\r)j)^ Let a be such a 
term and A^(a) the corresponding s-nucleus. By Ref. [24], Lemma 1, we have 



sign(A^(a)) = sign(a)(— 1 



|S+1 



Hence, all s-nucleus of G have the same sign if and only if all non-zero terms in the expansion 
of det(Z(i?, v) j) have the same sign. Further, there is an s-nucleus if and only if there is 
a non-zero term. It remains to remark that had we chosen any other set R of s reactions 
from TZ then there would be a zero column in Z(R,v)j {Isi~^o{Sj) = for all j) and hence 
det(Z(i?, v)j) = 0. This proves the equivalence between (ii) and (iii) when H2{i) = for all i. 

Denote the interaction network above by (A^cXq), Z{R, v)j by Zq and T{R)j by Fq. Assume 
that H2{j) is non-empty for some indices j. Any set R of s reactions such that T{R)j is non- 
singular consists of exactly one reaction among ^ Si and Si ^ for every i ^ J. We can 
assume that the reactions are ordered according to the species index i and let ii,...,ir be 
the indices for which outflow reactions are in R. Then det(F(i2)j) = (— 1)"^ det(Fo). For this 
set, if we change the sign of the columns ii, . . . ,1^ of det(Z(i?, v) j), the resulting matrix has 
a sign pattern smaller in the order ^ than that of Gj, i.e., of Zq. Hence det{Z{R,v)j) is 
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either zero or has the same sign as (— 1)"^ det(Zo). Hence, sign(det(Z(i?, v) j) det(r(i?) j)) = 
sign(det(Zo) det(ro)) or zero. 

Take a non-zero term in the polynomial expansion of det(Zo). It is of the form = 
=1= rii^j ■^o-s>5,,o-(i) with a a permutation on the set of indices satisfying i ^ J. Build a set R 
consisting of reactions — ?■ 5j if a{i) G Hi{i) and Si ^ (7{i) £ H2{i). Then, (— l)"^/! is a 
term in the polynomial expansion of det{Z (R, v) j) . This guarantees the existence of a set of 
reactions for which det(Z(i?, v) j) det(r(i2)j) 7^ 0. Consequently, M is injective over K.g{M,X) 
if and only if A/q is injective over ICg{Mo,Io) and the equivalence between (ii) and (iii) follows 
from the case above. □ 



Proof of Theorem [12^ Let a,b e W^, {X,w) £ ICg{M) andX = T(w). For each y ^ y' G 7^, 
define Sy^yi^i = and Vy^yi^i = for all i G ly^yi (that is, for all i such that Wy^y'^i = 0). Let 

My_,y> G M+ be such that My^yfa"^"^"'-' < 1 and My^y^b'^^ ''^^' < 1 for ah i G /+ , U /, 



y~^y y^y 



Then we can find Sy^yi^i G M+ and Vy^yi^i G M such that, 



V r ■ V / ■ 

, "^7-1^.-. =My^y'^?-^'" ^^^2^.^^,^^ =My^y,b:^-^'\ (A.8) 

"y^y',i + (^i '^y^y',i + "i 

Let d = {6y^y')y^y> (ZTZ, V = (v y^y>) y^yl ^-Jl , aud dcflne K = {ky^y')y^yl^Tl ky^yl = 

^y^y' l-^^v^v' where J is the cardinality of I~',,U/^ ,. Then S(v) = S(w) and = (K;,d,v) G 
K.h{M). Further, 



ly^y.a'^y^y' = ky^y, j| 



a/ 



i=i "y^y',t ^ "i 



and similarly for b. It follows that fx {a) = /K,w(a) and fxib) = /k,w(^)- This proves (ii). To 
prove (i) we follow the reverse procedure by choosing My^yi and Wy^y'^i to fulfill equation ( A.8 ). 

□ 
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